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PREFACE 

The application of artificial earth satellites in geodesy can be 
either dynamic or geometric » 

In the dynamic approach^ an evaluation of the perturbations in 
the orbital motion of a satellite can fbmish data from which the coef- 
ficients of the spherical harmonic expansion of the potential function 
can be determinedo Using this approach^ valuable information on the 
gravitational field of the earth can be obtained » 

In the geometric treatment^ the artificial earth satellite is 
used as a reference point or a triangulation station in the skyo In 
this way the positions of observation stations on the surface of the 
earth can be determinedo Some aspects of the geometric approach win 
be discussed in this paper* 

The purpose of this thesis is to serve sis a demonstration of the 
material studied during its preparation. I would like to emphasize that 
the reader should not construe the methods and procedures used herein 
as being of a standard type because of the following reasonsi 

1« The computed orbit from three observations at a 
given station at short intervals is assumed to be correct 
and accurate I 

2o It is assumed that the observed coordinates of the 
satellite taken fi*om Smithsonian Institution Astrophysical 
Observatory Special Report No, 66 and the computed orbital 
elements are referred to the same mean coordinate sirstemj 

3o The Gaussian Method used to determine the orbit is an 
approximate methodj 

iio It is assTomed that rally two approximations are necessary 
in the differential process of computing the geo^aphic latitude 



and longitude. 
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Mueller of the Department of Geodetic Science who has been ngr faculty 
adviser j to Mrs. Beatrice Miller of the Smithsonian Institution Astro« 
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Observatoryj to Dr. Nicholas T. Bobrovnikoff of the Department of 
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financial assistance and the opportunity to study Geodetic Science 
here at Ohio State; and to my fellow students whose thoughtful 
inquiries and interest in the subject served as a source of encourage- 
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1. INTRODUCTION 



Since the early days of history man has shovm a keen interest and 
an investigative concern about the earthy its characteristicsj and its 
dimensions o Since the third century B.Co when Pythagoras stated that 
the earth was of spherical shapes, many attempts have been made to 
measure its size and shape and to reduce the physical parameters and 
features on to maps and charts® Of the many existing earth sciences^ 
geodesy is the one that primarily contains this field® 

The inception of geodetic science can be traced back to Eratosthenes 
(circa 276 => 19it B«Co) who was the first man to measure the size of the 
earth® His method of measuring the length of a meridianal arc by timing 
the travel of a camel caravan between Alexandria and Syene (Aswan) and 
measuring the angle of incidence of the rgys of the sun shining into a 
well at these places was crude but effective® The degree of accuracy 
he obtained was remarkable when we consider the tools he had available® 
The science showed little development between this time and the 
seventeenth century® At this time the advent of the telescope, loga- 
rithmic tables and the method of triangulation (developed by the Dutch 
scientist Snellius in l6l^) took place and geodesy was reborn* 

The history of geodesy can be divided into periods based on the 
figure of the earth that was used as a reference during the respective 
periods® The studies and experiments of Pythagoras, Aristotle, and 
Eratosthenes brought the "flat”, or plane era to an end and introduced 
the spherical era® Newton and Huygens, in the eighteenth century, 
introduced the ellipsoidal period which fin ally phased into the 
geoidal and telluroifjal periods of this century® 
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With the launching, of the first artificial satellite and the begin- 
ning of the space age in this decade a new era in the field of geodesy 
is unfoldin^o By tislng. the geoid as a reference surface ire can determine 
the shape of the earth, but not its size e The orbiting artificial earth 
satellites iriU not replace the geoid or the telluroid but will provide 
the geodesist with a convenient reference station or **base-line in the 
sky" that can be used to measure both the size and shape of the eartho 
Since the launching, of our first artificial satellite (E^lorer 
I) on 31 January 1958, the field of space science has escpanded tremen- 
dously and with the successful results of the early and present projects 
the limits in this science are beyond imagination. At the present, 
the projects already xinderway and those planned for the immediate 
future suggest the foUcmlng. three areas of space activity; research. 



The applications include meteorologyj^ communications', and geodesy. 
No new principles, yet to be discovered are necessary 5 ,_no major tech- 
nological problems stand in our way. 

In the field of meteorology, the satellites can provide, through 
telemetry, information cm pressure, temperature, cloud coverage and 
stoim patterns and also assist in predicting such mahklllers as 
hurricanes, cyclones, etc. The day to day practical use of weather 
satellites is obvious. 

Xn communications, satellite contributicsis axe numerous: 
telephony, long-range radio ccmraunications, and international tele- 
vision. 





This brings us to the third application of space activity: 
geodesyo Geodetic satellites can serve both research and application 
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interests o As to the latter, navigation and surveying, come closest to 
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the immediate intei«sts of mankind. The geodesist can use the satellite 
in both a dynamic and a geometric approach. Dynamically, the perttir= 
bat ions in the orbits of the satellites caused by changes in the earth* s 
gravitational attraction can be mathematically computed and analyzed 
and thus provide valuable information to assist in determining the 
g:’avity field of the earth. By observing the changes in the orbits 
of many satellites over a long period of time, infozmation on regional 
gravity anomalies can be obtained. From these data ve can ccmpute 
the shape of our equipotential reference surface, the geoid. Further, 
observations of the nodal motion or the rate of change of the plane 
of the orbit of the satellite gives us information on the flattening 
of 1;he earth. Therefore, in the dynamic approach we can determine the 
shape of the earth. 

The size of the earth can be realized in the geranetric approach. 

In this process the satellite will serve as a reference point or a 
triangulation station in the sky. Thus we can connect triangulation 
nets or stations that are widely separated by oceans and continents 
or are located in inaccessible areas where the conventional methods of 
triangulation would be awkward and cumbersome. If so desired, the 
existing geodetic systems could be reduced to a common world geodetic 
network or datum. 

In this thesis a phase of the geometric use of the artificial 
earth satellite in geodesy will be discussed. In the geometric 
approach the geodesist can use the satellite in two basic ways. One 
method involves the observation of the satellite simultaneously ftcm 
two or more stations. The other method requires that the orbit of the 
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satellite be precisely known. Then knowing the orbital parameters or 
characteristics and predicting, their values for a desired time we 
observe the satellite from any munber of stations on the earth. 

The discussion of the problem will be treated in the following 



ways 

Firstp the orbit of the satellite Explorer I (1958 Alpha) will be 

computed ftom the right ascension and declination of the satellite as 

observed at three different times from the tracking station at Olifants- 

fonteinj, South Africa. The three observations have been made at very 

short intervals during, the same revolution of the satellite. The 

values of right ascension and declination are ones 

To tHV fPocH liso.o , 

from camera films of Baker°Nunn tracking cameras. The geographic 
coordinatesj latitude ( 0 ) and longitude ( ^ ) of the tracking station 
are known and its height (H) above sea level is given. 

Second^ the same satellite will be observed again from a station 
on the island of Curacao in the Netherland West Indies in the Caribbean. 
The geographic cooidinates of this station will be assumed to be unknown. 
Then using the computed orbital pstrameters and observing the satellite's 
right ascension and declination from this station^ the geographic co- 
ordinates of this unknown station will be found. 

Thirds the forward and back azimuths and the geodesic distance 
between the two stations will be computed. 

This connection between the two points is called a tie and in 
this instance $ an inter-continental tie. Ths distance between the 
stations is approximately 6^000 nautical miles, ^ho proeedu a re - dioouoood 
an tho following ■ ehapt e r s- oan - b e applied t o giy sat el lit e o b ee rv ed ftron 
■a ny . looation on tho ourfa e o of the ea rt h . 



2o DEFINITiaiS AJTO EIEMENTS OF THE ORBIT 



The motion of an artificial satellite in its orbit about the 
earth follows the same laws that apply to the movement of the planets 
in relation to the sun in the solar systemo 

Johannes Kepler (1^71 - I 63 O), a German astronomer and mathema- 
tician> derived the laws of planetary motion frcm the observations of 
Tycho Brahe (15U6 - I 60 I), Danish astronomer, with whom he worked as 
an assistanto 

Kepler’s Laws ares 

I The orbit of each planet is an ellipse with 
the sun at one of its foci. 

II The line (radius vector) joining the sun to 
each planet sweeps over equal areas of its 
ellipse in equal times. 

III The squares of the orbital periods of the 
planets are proportional to the cubes of their 
mean distances from the sun (foci). 

The relationship of an artificial satellite with respect to the 
earth is the same as that existing between the planets and the sun as 
stated in the above laws. 

In other words, we can replace the term planet vrith satellite 
and replace the sun with the earth and the laws are still s^jplicable. 
Now, since we know the basic form of the orbit is an ellipse we can 
easily describe the characteristics of the orbit. These orbital 
characteristics are known as parameters and can be considered to be 
constant quantities for a pairticulair satellite. The first two para- 
meters aire properties of the ellipse itself and tiiese are the semi- 
major axis, a, and the eccentricity, e. These two terms define any 
ellipse . 



5 



6 




Figure 1, 

Referring to Figure 1., the eccentricity is the ratio existing 
between the length of the semi-malor axis and the distance of the 
foci, c, from the center, o, or? 

e » c/a 

The relationship e^ - (a^ - b^)/a^, where b is the minor semi-axis, 
also exists. The eccentricity of a circle is zero and that of a 
parabola is equal to one. The eccentricity of an ellipse is less 
than one but greater than zero. The eccentricity of the orbit of 
a satellite is nearly circular, that is, slightly greater than zero. 
Two other parameters of the orbit are the angle of inclination, 
i, and the right ascension of the ascending node, Q . These two 
values define the plane of the orbit in space. Referring to Figure 2, 
the angle i is measured ftom the plane of the equator to the plane of 
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FIGURE 2 
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the orbit. It can have a value between 0° and 180°. The angle Q is 
measured eastward along the equatorial plane from the vernal equinox, 
to the ascending, node* It can have any value between 0° and 360°. 

The ascending, node is that point where the satellite crosses the equa- 
torial plane frcan south to north. 

Another parameter of the orbit is the argument of the perigee, 
w, or the angle measured firom the equatorial plane at the ascending 
node to the major semi-axis (the extension of which contains the 
perigee). The angle w is measured In the direction of the satellite's 
motion and gives us the orientation of the orbit in its plane. 

The sixth and last parameter is the time, Tq. Either the time 
of passage of the satellite through the ascending node or the time 
of passage through perigee may be chosen. The latter time is usually 
the one used in most computations* 

In summarizing! the a and e specify the size and shape of the 
orbit; i and Q specify the orientation of the orbital plane in space; 
w orients the orbit in its plane; T^ furnishes us with the value to 
locate the satellite in its orbital plane relative to the perigee. 

To position the satellite at aiy specific time we must resort to 
a system of angulnr values called anomalies, of which there are three; 
mean If, eccentric £, and true V« 
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Figvire 3. 



Refer to Figure 3» Angle QFN is the true ancaaaly V, which is the 

i 

angular distance of the satellite measured frcm the line of ^sldes 
along the plane of the orbit. Since the orbit of the satellite is 
an ellipse^ the radius vector FQ, or the distance of the satellite 
from the foci, will not be constant. Therefore, the angtilar velocity 
is not constant. This must be so if Kepler* s second law is to be 
fulfilled. If the angular velocity of the satellite was constant, 
the radius vector would sweep at an average rate and the eccentric 
anomaly would eq\:ial the mean anomaly. The mean anomaly is the theo- 
retical angxilar velocity of a satellite moving uniformly through 3^0° 
of arc per orbital period* 
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The Sobols and definitions of terms used in this paper are listed 
below ! 



a 



b 



c 



Semi-major axis of the orbit of the satellite 
(Chapters U, 5); 

Semi-major axis of the International Ellipsoid (earth) 
(Chapters 1, 2, 3) 

Semi-minor axis of the orbit of the satellite 
(Chapters U# 5)> 

Semi-minor axis of the ellipsoid (Chapters 2, 3) 
Distance from center of orbit to foci 



e 



f 



i 

k 

n 



P 

q 



r 



Eccentricity of orbital ellipse (Chapter 
Eccentricity of IntemationsJ. Ellipsoid 
(Chapters 1, 2, 3) 

Flattening of the ellipsoid 

Direction cosines (used together) 

Angle of inclination of the orbital plane 

Gaussian gravitational constant 

Ifean daily motion of satellite (velocity of the 
mean anomaly) in revolutions© 

Auxiliary parameter of orbital ellipse (replaces a) 

Geocentric distance of satellite at perigee 
(Chapter 5)s 

Time interval j (Ti + T 3 )k (Chapter k) 

Distance of satellite from observer 



s Geodesic distance 

t Time interval: (T 3 - T]^)k 

u Auxiliary parameter of orbit 

x,yj)Z * Geocentric space rectangular coordinates 



xai,ya,za Approximate space rectangular coordinates of 
the observer 



xoi,yo, 2 o Space rectangular coordinates of the observer 
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A,B,C,D,E,F,G Various terms used in Sodano's Fourth Method 
E Eccentric anomaly of satellite 

G Newton® s Universal gravitational constant 

H Elevation above sea level 

Ky Kl Gaussian determinants 

Kj (j = Constants for niitemational 

Ellipsoid used in Sodano's Fourth Method© 

L Difference in longitude between two points 

on the ellipsoid 

M Mean anomaly of satellite 

N Radius of curvature in the prime vertical 

R Geocentric distance of satellite 

Tj Time of observation (J = 1, 2j>3) 

To Time of Epoch 

V True anomaly of satellite 

X,TjZ Space rectangular coordinates of the satellite© 



CX Flattening of the ellipsoid (Ch^ter 6) 



Q( Geocentric Right Ascenion « the angular distal ce 

at the center of the earth measured eastward 
along the equator from the apparent vemal 
eqvdnox to the celestial meridian passing 
through the satellite© 

I 

Q( Apparent Right Ascension - measured at the 

position of the observer© 




Apparent right ascension obtained through 
the approximate coordinates of observer© 
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Forward Azimuth 
Back Azimuth 

Parametric latitrade (reduced latitude) 
Maximum latitude of the geodesic© 
Noicmal gravity© 
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Geocentric declination - the angular distance at 
the center of the earth measured northward or 
southward from the equator to the satellite, A 
satellite below the equator would have a south 
declination and the angle would have a minus value. 

Apparent declination - measured at the observer's 
position. 

Computed apparent declination obtained through 
the approximate coordinates of the observer. 

Spherical distance 

True spherical distance 

Greenwich hotir angle of the apparent vernal equinox , 

Local hour angle of the apparent vernal eqviinox. 

Geographic longitude (positive westward and 
negative eastward) 

Geocentric distance of observer 

The angular motion of the satellite during the 
interval between observations. Usually expressed 
in radians. 

Geographic latitude (north positive and south 
negative ) 

Geocentric latitude 

Argument of the perigee. The angle measured 
eastward in the orbital plane from the line 
of nodes to the line of apsides. 

G>aussian determinant 

Vernal Equinox ” the point of intersection of 
the path of the sun with the equator at which 
the sun is passing from south latitude to north 
latitude « 

Angle at satellite's position measured between 
the line to observer's position and the line to 
the earth's center (parallactic angle at middle 
observation) . 

Angle at observer's position measured between 
the line to the earth's center and the line to 
position of the satellite at the time of the 
middle observation. 
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Right ascension of the ascending node • 
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3. DETERMINATION OF THE MASS OF THE EARTH 
AND C5ADSS* GRAVITATIONAL CONSTANT 



The existing relationship between the earth and the satellite 
can bo treated as a modified two body problem in which we assvane that 
the two bodies are spheres and that their mass distribution is 
homogeneaiSo Under these circumstances, we can say that they will 
follow the laws of celestial mechanics and gravitation. That is, they 
will attract each other with a force that is proportional to the 
product of their masses and which varies inversely as the square of 
the distances between their centers. Since the mass of the satellite 
canbe considered to be negligible compared to that of the earth and 
because of its proximity, it cante stated that the gravitational attrac- 
tion of the earth is the only attracting force influencing the orbit. 

The force acting between the two bodies is: 



where r is the distance between the two bodies, me is the mass of the 
earth, ms the mass of the satellite, and k equals the Gaussian gravita- 
tional constant. The numerical value of the gravitational constant 
depends on the units of mass, time, and distance chosen. ^2, P.^?] 
Kepler's third law: 



is a simple relation between mass, time, and distance that may 
conveniently be used in practical application of elliptic motion to 
find any one of these three quantities when the other two are known. 
The units of mass, time, and distance can be arbitrarily chosen and k 
can either be computed or derived from observations. In applying the 
xmits to determine k for use in observation of planetary orbits the 
following units are used: the sun's mass as the mit of mass, the 



F s lc2 me mg 




k^ (me + mg) = n^a^ 



f j j 

m 
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ephemeris dapr as the mit of time, and the semi-major axis of the earth's 
orbit (which approximates very closely, one astronomical unit) as the 
unit of distance o In the case of the artificial earth satellite a 
revision of these units must be made in order to find a realistic value 
of the Gaussian constant that will be applicable. Bear in mind that 
we are theorizing that the earth is the only body influencing the motion 
of the satellite. This is a sound theory but the ether assumptions we 
made - that the earth is a sphere and that the mass distribution is 
homo^neous” are not absolutely correct. In fact, if the earth was 
spherical and the mass distribution was constant throughout, then the 
orbit of the satellite would be nearly perfect and very simple to 
predict. The plane of the orbit would be fixed, the parameters of the 
elliptical orbit woiild remain constant, and the center of the earth 
would act as one of the foci. 

However, the earth is an ellipsoid of revolution with a bulge of 
mass at the equator and a non-homogene ous distribution of mass through- 
out. These characteristics plus the effect of the earth's atmosphere 
and other perturbations' have a marked influence on the behavior of the 
satellite in its orbital travels. 

The effect of the attraction of the equatorial mass of the earth 
is to pull the satellite toward the equator, thus decreasing the inclin- 
ation. Behaving like a gyroscope the actual effect is to cause the 
plane of the orbit to process westward for an eastward launch about the 
focus. This motion causes the perigee to advance to different points 
along the orbit. The perigee moves in the direction of the satellite 
motion when the inclination is less than 63.°5 snd moves inttie opposite 
direction when the inclination is greater than 63.°5» Consequently, 
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the elliptic orbit rotates about the focus in the orbital plane. Reduc- 
tions of the periods of this peidgee motion from the observations provide 
data to determine the shape of the earth* 11, P»U 



The non-homogeneity of the earth’s mass distribution causes per- 
turbations in the orbit of the satellite. By observing a nvunber of 
satellites for long durations of time we can compute the deviations 
from the normal gravity and thus obtain values for the regional gracvity 
anomalies* 

Proceeding vriththe computation of the Gaussian constant, k, the 
following units will be used: the eairth’s maiss as the unit of mass, 

the ephemeris day (86,iiOO seconds) ais the unit of time, and the semi- 
major axis of the earth as the unit of distance* Compare these with 
the Tonits mentioned earlier in this chapter that are used in conjunction 
with planetary orbits* 

Note that aU computations are based on the International Ellipsoid 
(Hayford, 1910) which has an equatorial radius of 6,378,388.000 meters 
and a flattening of l/297»0. Here again rises an arbitrary value* 

There are many reference ellipsoids from which to choose but the varia- 
tion in size fi*om one to another is small* The International Ellipsoid 
has been adopted by many countries and it is the one that is used by 
Smithsonian Astrophysical Observatory in their computations, so I 
decided to use this as the reference body. 

The intemational gravity formula adopted by the International 
Union oj 



where 




- gravitational flattening, € ® theoretical constant* The 
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coefficientts of the formula are: 

7 = 978,01:9(1 + 0.0052881: sln20 - O.OOOOG59 ain^Z(:j))c Vse<?^ 

The gravity value at the equator is the resultant of the gravitational 
attraction and the centrifugal force due to the earth's rotation. 
Therefore, 

Te ” ir^a * 978,01:90 cm/sec^ 

where Te “ attraction of the earth at the equator and rr^a - centrifugal 
force, 

27T 

If rr^a ~ (-r^)2rjwhere w is the angular velocity of the earth, 

T is the sidereal day in seconds^ and r the radius of the ellipsoid then: 

= 3.38051:5 cm/sec^ 

Then 

Te - 3.3805 = 978.01:90 

Te = 978,01:90 + 3.3805 = 981.U295 gal 



where a gal (from Galileo) equals a unit of acceleration equal to 
1 cm/sec^. 

The mean density of the earth is computed firom the formula ^5] s 

r. = 2TTe (arctan e' - ) 

where P© = Te/ a - attraction of earth at the equator divided by the semi- 
major axis, G s Newton’s xmiversal constant of gravitation which equals 

6,673 i 0.003 X 10”® cgs, e - the secondary eccentricity, and (Xm = 
the mean density of the earth. To solve the equation, two of the three 
quantities Te, a, or 0^ m must be knovm. In our example, the Te has 
been determined aid the value of the semi-major axis of the Interna- 
tional Ellipsoid is used. Tte computed mean density is (see Table 1): 
CTj. s 5.5130OU37 gm/om^ 

The volume of the International Ellipsoid is computed from the 



formulas 



t 
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V = 1+/3 7 T = 1.083320781 X lO^"^ cnp^ 

where b - 6,356,9U*9U62 meters, the semi-ndnor axis of the ellipsoid. 
Using, these values, the computed mass of the earth is: 

/v; r \Z(/ =(1.083320781 x lO^”^) (5.^1300^37 )= 5.97235220 x lO^^gm. 



To obtain the proper tinit for k in terms of mass of the earth for 
vinit mass, time in days, and distance in terms of the earth’s semi- 
major axis the following ccmversion is made: 



sec* 



■ cm- 



.3 = 



cm*^ 



m 



^ r, (86,U00)2 (5.9723522 x 10^^) = k 

(6,3-/8,3tiaoop 



and k = 107.0731376 and k^ = llU 61 i .6568 
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TABIE 1 



CCMPUTATIQN OF VOLUME, DENSITY, AND MASS OF THE EARTH; CCMPUTATICN OF k 



DIMENSIONS OF THE INTERNATICMAL ELLIPSOID 



CC “ 1/297 

a - 6.378,388.0000 meters 
b - 6,3^6, 911*9U62 meters 
e = 0.0819918898 

ef = 0.0067226700 

e . = 0.0822688896 

• 0.0067681702 

= 0.0005568098 



1 + ^2 



1.0067681702 = 

“ 700 S 555 W ■ 



1808. 10066 



ARCTAN e' = e' - ^ “ £1 + • • • 

... = O.O8208I1O365 

/ 

^ = 0.0817158230 

1 -f 



r ( 2 X 3.1^159265353^ 6,378,38800. = 3.3805U56U00 

T y 6 ,uoo 

M • h/3TT a^b = 1.083320781 x 10^7 cm^ 

= Te 



27 rG 



1+e' 



7 T 



5.51300U37 gpi/cm^ 



e ^3 



arctan-e' — . 



1 + e- 



n 



(a) 



M = VCfsr 5.972352200 X lo27 gm 



k = ^ .t cm^ - (J 
se^ 






gm sec‘ 



gm = G (86.1tOO)2(5,972352 x 10^7 

(6,37&,38800p 



107.0731376 



* 






U. DETERMINATION OF THE CRBIT 



Since the or bit of the satellite is defined by six parameters then 
six independent quantities must be given by the observatims in oi^er 
that the parameters may be determined. A single observation gives two 
quantities, the angular coordinates of the body. Therefore, three 



The apparent positions of the observed satellite are obtained 
by measuring its angular distances and directions from a background 
of stars on the film of the Baker-4Iunn camera. The star9 on the film 
can be identified from catalogues containing their right ascension and 
declination. With these coordinates known, then the right ascension 
and declination of the satellite can be cranputed. Catalogues containing 
the observed values of time, right ascension, and declination reduced 
from the camera film are published about six months after the observa- 
tions were made. These reduced values are quite accurate. However, a 
more precise reduction of the films is made by the Smithsonian Institu- 
tion Astrophysical Observatory Photoreductim Division. These results 
are published and made available about 18 months after the observations 
are taken* 

In this paper the con^mtation of the orbit of Explorer I (1958 Alpha) 
satellite will be reduced from observations of right ascension and 
declination from the Baker-Nunn Camera Tracking Station at Olifants- 
fontein. South Africa. Date of the observations was 18 December I960* 

The observed data I?, p.lO sr 



complete observations are sufficient to define the orbit. 




Time (UT) 



(X 




22^ 29 10® .m 

22^ 31“ 0itS*li2 
22 “ 32 ®“ 03®. 10 




-U7° 22'00« 
06’00" 
- 1 * 9 ° 2 ^' 00 " 



These data are not those precisely reduced by the Photo-reduction 



I. 
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Division but the precision of the observations is as follows: 7> P«3] 



Standard Error in timing ( ) 

Standard Error in direction ( ) 



. 00 ^® 

2' .7 



gt ^ .02® 

o; = 3'. 5 



'T> = 

The geographical coordinates and hei^t of the tracking station are 



[8, p.7] 



0 = 25° 57” 3V'. 70 S 

A = 28 lU' 5l'M0 E 

f-f = 1514; meters 

There are two classical methods available for determining the orbit 
one developed by Laplace in 1780 and the other by Gauss in 1801, Both 
methods were devised for use in determining the orbits of comets and 
planets about the sun but are also applicable for tise with artificial 
earth satellites. Since the development of these two methods many 
papers have been written on the theory of orbital determination but 
very little that is theoretically important or radically different 
has been added to the methods of these two men. The Laplacian method 
is not as complicated ais the one devised by Gauss but it has two weak- 
nesses. It is inaccurate when the Intervals between observations are 
short and it assumes that the observer is in orbit about the focus of 
the ellipse. The Gaussian method, even though it is an approximate 
method was chosen to determine the orbit. There are other more modem 
solutions available but these are not as accurate. 

The first step in the calculations is the conversion of the 
Universal time of the observations to sidereal time. This is required 
in order to introduce "the vernal equinox into the coordinate system. 

The sidereal time at any instant is the hour an^e of the vernal 
equinox. 

The origin of the coordinate system is at the center of mass of 
the ellipsoid. The x -axis points toward the vernal equinox, the y-axis 






m 



* 



I 








I 



» 



. I 

I • • 




f t - 
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is directed toward a point 90® eastward from the vernal equinox, and 
the z-axis coincides with the north pole of rotation. The x-y plane 
lies in the plane of the equator. 

The rectangular space coordinates of the observer are computed 
from the given latitude and longitude of the tracking station. Then, 
by changing the sights of these values the coordinates of the center of 
the earth are obtained. The observer’s position is computed as follows: 

Xj = (N + H) cos 0 cos Q (j = 1, 2, 3) 

yj = (N + H) cos 0 sin Q ij 

zj = (N + H) sin 0 (1 - e^) 

The position of the observer, the earth* s center, and eventually the 
position of the satellite are all converted to rec(tang\ilar space 
coordinates to provide a more convenient and a common system of 
coordinates. 

The next step is to compute the direction cosines from the 
observer's position. These are functions of the observed (appsirent) 



right ascension and declination and are obtained from: 
f ^ = cos 0 -5 cos (X . 



U = 1, 2, 3) 



gj = cos (5 j sin CX 

hj = sin 



,3 

3 



The orientation of these axes: f-axis points toward the vernal equinox, 
g-axis perpendicular to f-axis, and the' h- axis parallel to the axis of 
rotation of the earth and i)erpendicular to the f- and g-axes. The direc- 
tion cosines give us a line in qpace from the observer'^ posi’^ion to the 
satellite. The computations of these three steps for the three observa- 
tions are contained in Table 2, 



TABLE 2 
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CCMFUTATICN OF SIDEREAL TIME, COORDINATES OF EARTH GEMTER, DIRECTION COSINES 



Date 

Sidereal Tiiae At O" UT 
UT Time Interval 
UT To Sidereal Time 
Greenvrich Sidereal Time 



18 December I960 
^H6-25.U19 

22-29-10. lliO 
3-U1.633 
■5U-19-If35'2 



X (Add East, Subtract West) +lr^2-^9.it07 
Local Sidereal Tine (19 Dec) 

0 ^ 93°0ii«08".985 

sin (/ T .998^6561* 

-.0535U139 

5^-07 “^-o 6® 
76°ii6'30« 
.973li7917 
.22877565 

J*7°22»00» S 
-.73570317 
.67730UIO 



1 - e^sin^f 
(l-e^sjm^C 
a/d-e^Bin 



ii i 

<p )I:N 



18 December i960 
5-i^6-25.Ul9 
22-3I-OU.U20 
3-ltl.9lt6 

^-21-11.785 

+ Ir52-59.1i07 

“StHiiSlWr 

930321 U7". 880 
.99808U78 
-.06186091 



7^.28^-514® 

H 2 ° 13 ‘ 30 '' 

.92570563 

-.3782l4i47l4 

5 l° 06 » 00 » S 
-. 7782 U 315 
.62796306 

25 ° 57 ' 3 U ''.7 S 
-.U 3773790 

.I916H4I47 

.89910263 

6,378,388.0000 ^ 

.00672267 ^ 

.99327733 

.00128816 < 

.998711814 

.99935571 

6,382,500.18 meters 



18 Dec. i960 
5-I46-25.I4I9 
22-32-03.100 
3-I42.IO7 
28-22-r(576’55 
+ Ir52-59.l407 

93°l47'30«.l495 

.9978109I4 

-.06613118 

8l\-33^.00® 

128°15*00» 

.78531693 

=.61909395 

l49°25'00» S 

-.759I46058 

.65055333 



H 

I^a 

N/a 

N/a + H/a - N+H 





I5l4l4 

.OOO 2 I 42 O 7 

I.OOO 6 I 4 I 47 I 

1.00088678 


meters 


^ 








> 




-(N + H) cos^^ cos = X 

-(N + H)cos (n sin * y 

-(N + H)sin ^ (l-e2) =z 


+.OI 48 I 8 I 89 
-.89860916 
+ .I 4351807 O 


+ .05566863 
-. 898 I 76 I 43 
+ .I 435 I 8 O 7 O 


+ .0595111414 
-.89793000 
+ .I 43518070 


cos 0 , cos Of' - f 
cosO sin Of' = g 
sin^' = h 


+ .I 5 I 495069 
+ . 6593 I 1 IU 3 
-.73570317 


-.23752372 
+ . 58 I 3 O 89 I 4 
-. 7782 I 43 I 5 


-.I 40275363 

+. 5108905 I 4 

-.759146058 



( 



2U 



The next step is to determine the distance of the satellite from the 
observer's position and from the center of the earth. This is done by 
solving the triangle formed at the time of the middle observation between 
the center of the earth, the observer, and the satellite. To solve the 
triangle the satellite is "stopped" in its orbit at the time of the 
second or middle observation, (T2)* 



The geometric situation is depicted in Figure Since SOE form a 
triangle at T2 the quantities r 2 and R2 must satisfy the equation 



This equation treats the geometric situation. In the Gaussian method 
a formula relating the two unknown quantities r and R applies dynamically 




s 
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where K, A , and are determinants involving the coordinates of the 
earth’s center and the direction cosines from the observer's position. 
These determinants and their values are: 



A 



K 



4 



f 1 f 2 f3 
61 62 63 
hi h^ h 3 


S 


- O.OO 20 U 368 


' ^1 
^1 


XI + X 3 

7l + 73 
zi + Z3 


- 2 X 2 

- 2 y 2 

- 2z2 


f3' 

63 

h3 


= + 0.00018897 


fl 

61 


XU- X 3 

7l+ 73 

21 + Z 3 




r 1 


= -O.U57l61t31 



k'i= - 

The solutions of these determinants are shown in Table 3o 

The two \mkncwmsr2 and R2 are expressible in terms of a single 
unkncfwn, the angle $ , by means of the equations ^U,p.235] : 

Pg cos 2 = X2 f2 + y2 62 + ^2 ^2 ( 0 ^ = 7T ) 

N sin m - P2 sin 2 ” .U8507360 

f\| cos m = P2 cos ^ o - JL- « -.78l55lii6 

\iy A 

The computed values for I 2> N, and m are: 

^1/2 = 150° 58' 12".06U 
N = Z 0.91981^731 
m = llt8° 10' 26«.395 

The solutions for these quantities are contained in Table 
N, m, and M are intermediate parameters used in the solution for the 
angle $ , where [it, p*J12] : 

li= + 0.163U.533 

The sign of N is chosen ao that M shall be positive. 

The equation for the solution of the angle ^ is [I4, p.213] : 
sin ^ = M sin ( (^ + m) 

which is a transcendental equation. The quantities M and m are known and 
U is positive. One solution of the problem results in the position of 



( 
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SOLUTICN OF DETERMINANTS 



A 



f, g2 h3 = - .068 U07 82 

^ g3 hi - .089 276 57 

£3 Si “ .206 66 U 15 

f3 g2 hi - .172 21*5 98 

fl S3 h2 * - .061 607 9h 

f2 gi h3 - .118 938 5U 



K fl (n + 73“2y2)h3 

(X1+X3 - ^2)s31^ 
^Si(zi+ Z 3 - 2 Z 2 ) 

5 (71+73-272) hi 
fi S3 (zi+Z3-2z 2) 
(xn- X 3 - 2 X 2 ) gih 3 



.000 021 92 
.001 369 62 
0 

-.000 055 20 
0 

.001 82li 68 



Kl 



fi (yi+y3> J3 

(XI+X 3 ) g 3 hi 

^3 SI tzi + Z3) 
£3 (71 + 73) hi 
fl S 3 (Z 1 +Z 3 ) 

(X1+X3) gi h3 



.211 


1*11* 


83 


-.01*0 


1*78 


02 


-.231 


126 


31* 


-.532 


327 


35 


.068 


900 


28 


-.053 


926 


77 



1 



27 



TABIE U 



SOLUTICW OF%, N, m. 



xj f 2 

7 p So - 

z| h| = 

X 2 f 2 + 72 g 2 -*- *2 ^2 = 



-.013 222 62 
-.522 117 99 
-.338 676 UO 

-.b7ii 017 01 




cos 






2 





.003 099 00 
.806 720 90 
.189 382 2 h 

-imiunn 

.999 600 99 



X2^2 + 72g2 + Z2^/P2= 



Bin 






P2 sin ^2 K- 
P 2 'j;/ 2 - 

sin2 m 
n 2 cos 2 m 
N 2 (sin 2 ni+ cos^m) 

N 



- .871 365 89 
150 ° 58 » 12 «. 061 i 

.W5 267 23 
.U85 073 60 

-.781 551 U6 

.235 296 UO 
.610 822 68 
.8u6 119 oB 
.919 6 h 7 31 



P2 sin 
f/p ®°s 




.527 3U1 1;3 
-.8U9 653 hi 

lii8° 10»26'«.395 



A 



JOl 



sin^M/ 

N 

t 2 

U A p s±n^4/ 
h APi Nsi^^ 

./b 2 Ki 



Solution for t: 



t = 



-.002 0 U 3 68 

.998 803 lt5 

.nU 272 81 
.919 8Ii7 31 
.011 485 88 
-.457 164 31 

-.000 858 24 

-.005 250 93 
.163 445 33 



k (Tj - Ti) 



P 



2 



cos 






2 



N sin M 
N cos U 



sin m 
cos m 



M 



T3 = 22h 3^ 03s. 10 

Ti = 22^ 29 ^ IOS.14 

^ - Ti - ^ 52 ^^. 96 

I3 - Ti = .002 001 85 days 

k = 107. 073 1376 

k(T3 - Ti) = 0.214 344 36 = 2 t 

t = 0.107 172 18 



t 

f 
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the observer and is not applicable. That is, ifr = 0, P = R, The poxiper 
solution for the angle must satisfy the inequality [ii, p»2l3] : 

(f) - 7T ^ ^ 

In our problem, ^ equals l^O*^ 58'12’',06U so ^ must be less than 
29°02'. The solution maybe obtained graphically be plotting the curves 

[U, p.213] : 

yi » sin^^ T 
72 = M sin ( ® + m) 

Once the approximate value is found from the intersection of the two 
curves then the exact solution may be found by numerical trials or by a 
differential formula. The solution for the angle is lengthy and time - 
consuming. The approximate value from the intersection of the curves on 
the graph was found to be approximately 23°. The exact value, solved for 
by numerical trials was: 

$ - 230 1^‘ I49”. 16 



Hith this value knovm we can now solve for the two unknown sides of the 
triangle SCE in Figure li, which will yield the distance of the satellite 
from the observer’s position and the distance of the satellite from the 
earth’s center for the second observation. The formulae used and the 
results are [^U, p. 25^j : 

Rj = Pp sinTz . 1.228 150 03 

$ I 

r, . Pp 5ln( 'i'-*- )- 0.25U 280 87 
— 



The observer - satellite distances for the first and third 
observations are given by ^ii, p. 2^5]: 



nf3 

glg3 



ri = 



xif3 

yi&3 



■fei) 



^2^3 

y2e3j 



(1^2) r 



,^3% 

[7383 






f2f3 

8283 



I 
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flf3 

g-lg3 



^3 




2173 




h^2 



yrherB 







2.93U 11079 



1.9U3 20695 



O.51U61323 



Cl>- 3 ) = 1 = 1.50993233 

tCtJ I q . 

2 ‘ 2t + IjR^ 

where 

t = = 0.10717218 

q = K I -0.03UU5185 



The time of the second obsesrvation, T2, is the origin of time. Upon 
substituting these values into the above formulae and solving the 
determinants, the observer - satellite distances for the first and 
third observations are found to be; 



ri = 0.23531756 
r3 = 0.27239982 



The detailed solutions for Ra> T2, and q are listed in Table 5. 

The next step is to determine the geocentric space coordinates of 
the satellite for the time of each observation as follows p.256j : 
r 0.23531756 O.25U28087 0.27239982 

X = r f - X 



r = r g - y 



-0.01171927 -0,11606637 -0.1692211*5 

+ 1.05376378 + 1.01*599217 + 1.0370961*9 



I 



TABLE 5 



DETERMINATICN OF R2, V2, q, ri, r^. 





Pz 


sin 


w 

T ^ 2 




O.U85 073 60 




Pz 


sin( 


% 2/sii.'P 


3 


0.39li 962 82 




sin 


cx> 


1.228 150 03 = 




'V, 






s 


150 ° 58 '12 ". 06 




( % + 
sto J 




8 


_^Q.3.5.Mi9".a6 
17U^ 1U»01".22 

0.100147158 




. Po 


sin 


{'i'f+4) ) 


m 


0.100li311i9 




P 2 sin( ^ 


2 + 1 |) 


1 — 


.251i28087 = 


T 3 


• 21 


32 


!“^ 03 ® .10 




T2 * 22j^ 31 “ OUS.1^2 


T2 




3r^ 0 I 1 S.I 42 




Ti = 22 “ 29® 10®.1U 


T 3 - T 2 


= (+) 


58®.68 


T2 


-T]_ = (-) 1® 5U®.28 



T2 - Ti - (-) 

- Tr, a (+) 58® .68 

Ti + ^ a (-) — 5505 
Ti + To (days) = -0.0006U352 

k = 107.0731376 

k (T1+T3) = -0.06890371 

k (Ti+T3)/2 = - 0 . 03 UU 5185 = q 



t 0.10717 218 

2t O.2IU3UU36 

t2 0.011148588 

Ro 1.22815003 

R2I 1.852U8316 

2R2^ 3.70U96632 

1 - q/t + t2/2E23 = 




l4R2^ 7.U099326U 

q -0.03Ui5l85 

q/2t -0.16073131 

q/t -O.32IU6262 

t 2 / 2 R 2 ^ 0.00310013 

t2/l4R23 0.00155007 



1.32U56275 1+^t + t2/2R23 * 0.68163751 

X|^| j ■» 2.93lill079 



1 -I + ^ ° -fel] ° 1*9W20695 



1.50993223 



TABLE 5 (Continued) 
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S3 

n £3 

^is3 - nf3 




0.02li6l^7 

0.36m8l0 

-■53575C2H5 



X2E3 

y2^3 X2f3 

^2^3 ~ ^ 2^3 “ y2&3 



0.028U*058 

0.36l7l;382 

- O.3333032I4 




- 0.97791863 



^3|3 

X3g3 - y3f3 





X3f3 

y3t3 



0.0301^0383 

0,36l6l4li^7 
- O.3312I167I4 



- O.6U36693I 



^2^ 

22^3 

^263- 82^3 



f2f3 

82^3 



^2 




^“2^3 

§2^3 



- O.I213U862 

- O.23UI2U29 
0.11277^67 



0.08U1U060 



^183 

Slf3 

flga*” 81^3 



fifa' 

glg3 



0.0791628U 

0.26555215 

0.3¥i7llt^9 



^1^31 

8183 J 



0.081117U9 



= .081117U9/.3Ui71i*99 = 0.23331756 



^■lyi 

Sl^i 

f m - ei^i 



tixi 

si^ 




fixi 

81X1 



- 0.1392ii011 

0.03176832 

o.inooBUa 



- 0.08800320 



glX2 flX2 

f 1 X 2 - 81 X 2 = gl32 




^■1^2 

81X2 



0.13917306 

O.O3670U63 

O.I75B7755 



- 0.26556339 



V 



TABLE 5 (Continued) 



^U3 

glX3 

^iy3 “ 81^3 



flX3 

8173 



- oa3913i;87 
0.03923836 

- O.17B37323 



fl|2 

21^2 

^182 “ &1^2 



flf2 

8182 



O.O9OO7I422 
- 0.1^660923 
0.2li6683U^ 



T2 




flf2 

8182 



0.09U7133U 



^1^3 

8163 



0.09390030 



- .09390030/*3UU711i99 = 0.27239982 
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2 = r h - z - 0.60830U57 - 0.63307305 - 0.61*205763 

1^ ^(x^+^+2^) 1.216791*93 1.22815003 1.2311*391*0 



There is a check for acciiracy of the computations Txp to this point 
in this last step. Note the agreement in the value for R 2 computed from 
the space coordinates and that determined fi*om the formula, (P 2 sin '^2)/ 
sin ^ , in Table 5* 

The inclination of the orbital plane and the intersection of the 
orbital plane with the plane of the equator (XT plane) can now be solved 
using the geocentric space coordinates in a ratio of determinants. This 
solution will produce the first two parameters of the orbit that we are 



seeking: i and Q [u. 


p.256 


1 = 








A. n. c = 


Yi Zi 


• 


'ZiXi 


• 

• 


XI Yi 




Y2 Z 2 




Z 2 X 2 




X2 Y2 



See Table 6 for the evaluation of the ratio of these determinants and 
the detailed solutions for- i and Q • 

Then, 



cos 1 » C 





sin Q sin i = 


A 




^a2 + b2 + 


cos ^ sin i - 


-B 



+ b2 + 



i « 32° 3l*»20«.2U5 



Q = 206° 00* 27". 603 



If the inclination was greater than 90° the sign of -the previous 
two equations would be reversed. The values i and ^ uniquely determine 
the position of the orbital plane in relation to the equatorial plane 
of the earth |^U, p«ll*6j* 

The next step is to find the auxiliary parameter u for each observa- 
tion. Solution of this value is necessary in order to find the "longitude" 



TABLE 6 



EVALUATION OF DETEraHNANTS As B: C j SOLUTION FOR i 



Ti 22 

I1Z2 - I2Z1 



Z1X2 

Z2X1 

Z1X2 - Z2X1 = 



X1T2 

X2Y1 

X1Y2 - X2Y1 « 



YlZl 

Y2Z2 

L^ 2 ^ 2 J 



Xlll 

X2Y2 






a 2 + b 2 + c 2 




a2 + b2 + C 



/V 



sin 1 



a 2 + b 2 + 

A / 



V 



^a2+ b 2+ c2/sin i = 



bA /a 2 +b 2 + c 2 



_E_ 



-^a2 + b 2 + C^siji i" 

Q 



-0.667109li5 
-O.63628I82 
-0.03082763 = A 

0.07060370 

0.007U1915 

“o7o55IHII55 = b 

-0.01225826 
-0.12230656 
o.iloou82l] ® c 



0.13058811 

0.8U271286 



32°36’20".2U5 

0.53836329 

-0.23606766 



- cos i 



-0.63869138 = sm 

0.68386612 

-0.89873535 = cos 

206°00«27".603 



inQ 



Q 



3U 

4 ND n . 
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of the satellite from the position of nodal passage. The longitude of 



the satellite from the node is called the argument of the latitude. 




Figure 5« 

Refer to Figure 5® The formulae that define u^, U2, u^ are |^ii,p.238| 
cos Uj = cos bj cos (Ij - 0) (j«l, 2,3) 

sin Uj cosi = cos bj sin (Ij - Q ) 



sin Uj 



sin b 



y/si^i 



To determine the parameters a, e, and w^we must solve for the ratios 
of the triangles formed at the time of each observation. With the geo- 
centric coordinates and the range of the satellite from the earth’s 
center kncwn the above equations may be substituted irith the following 



expressions 



k , P.253,257] 



sin i sin Uj 



cos i sin uj 



cos Uj 



I 

- Y 






cos 



— 0 
“ sin ^ 



j * (1, 2, 3) 
- fL sinQ 



+ cos 

RJ 



Then for abbreviation? 

a U3 - U2 = 2°.iUa97583 = 0,0U262050 radians 

= U3 - ui » 7°.U5516 o83 » O.I 30II7IO radians 

(Ij - U2 - ui s ^o, 0131S500 = 0.087U9658 radians 

Then the ellipse parameter p is determined and used in place of the semi- 
major axis, a. 



i 



i 

I 
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The equation that defines p is p,257j ; 

k ^ Cf^ + R-)^ ^ R~^ 02 ^ ^ ^ “• 3 ) 

60'i6'3 6cT3 6 (f^ 

Using this formula the value for p was found to be: 

p = 0.d30^$3$h 

With the values of p and u now known the eccentricity, e, and the 

argument of the perigee, w, can be found from |^U, p,257] s 

e sin (ill “W) = R3(p~Ri) cos (^2 - Ri (p *Rt) 

Rl R3 sin (^ 2 

e cos (ut-w) = p - Rl 

Then the semi-major axis of the orbital ellipse is defined by 
® • '{'1 - « 2 ) ' 

The computed values for these three orbital elements are: 

e = 0.32815828 
a = 0.93078812 

w = 82°51i'ii2».337 

The detailed computation of these parameters and the auxiliary para- 
meters U2_, U2, u3, aid p are listed in Tables 7> 8, and 9. 

The next and final orbital parameter to be computed is the time of 
epoch. To (time of perigee passage). This value is computed from Kepler’s 
equation: 

M - E - e sin E 

after the true and eccentric anomalies are determined for the three 

observations. Values for To ares 

Observation #1 Observation #2 Observation #3 

December l8«9Ui5U9i;2 I8.91ii53l82 I8.91ii55056 

Detailed computations of Tq are listed in Table 10* 

To verify the results of the computations, the orbital parameters 
were determined again using data published in Smithsonian Institution 



37 



Astrophysical %servatory Special Report #78 dated 25 October 1961, This 

publication contains both the mean elements and, the smooth elem«its of 

Explorer I (1958 Alpha) satellite based on 82 observations made fi^jm the 

various Baker-Nunn Tracking Stations for the period l6 December through 

29 December I960, The published smooth elements ares 

To “ 37291.0 Modified Julian Days (Smithsonian Days) 

w = (lii0°,6 ±2)+7°.275 (t-To) + 0°,000ii27(t-To)^+0°.22li cos w 

Q = (1860.278 ± 9)-U°.896(t-To)“0°.000298(tuTo)2+0O,00887 cos w 

i = (33°.210 ±3)“0°.00055U(t-To)+0°.50xl0“^(t-T^)2 -0°,OOU73 sin w 

e = (.09971 ±8) -0.000202(t~To) + 0.30 xlO“'7(t-To) -t-O.OOOii55l9 sin w 

M = (,9liOl6 ±8) (13.1i35l3U ±5) (t-To)+ (0.0002839 ±6) (t-To)^ 

-(.11 ±1) xl0-5(t-Tc,)3 -(.000631 xlO- 3 ) cos w 

Standard error of one observation; m = i5’’®U3* 



The angles are expressed in degrees, e in terms of the earth's semi-major 
axis (6,378,388 m,), and M in revo3.utions. The Smithsonian Day for 
December 0, I960 is 37268.0. Using the time of epoch, Tq “ December 
I8,911i53l82 from the middle observation of our calculations and redesigna- 



the term (t - Tq) will have the following value; 



ting it ais t 

Date 

December 0, i960 
Date of Observation 
Date of Observation (t) 
Epoch (To) 

(t -To) 

Detailed computations of the elements 



Smithsonian Day 

37268.0 

+ I8.91li53l82 
37286.91h53ia2 

37291.0 

-l4..0b5U681d days 
are contained in Table 11. The 



computed and published parameters are; 



Published 



Computed 



i 33°13'00.1i02 

Q 206°16'U3".067 

e 0.10011028 

1.17098895 
HO°57'53”.75U 



32o3li'20.2U5 

206°00'27"*6o3 

0.32815828 

0.93078812 

82°51i*ii2".337 



a 

w 



TABIE 7 
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CCMHJTATICW OF THE AUXILIARY PARAMETER, U 





Observation 1 


Observation 2 


Observation 3 


Z^/R. 

Sin 1 

(Zjy^j)/sin i - sin u 


-O.U9992366 

0.53836329 

-0.92859909 


-0.515U6882 

0.53836329 

-0.957U739U 


-0,52138792 

0,53836329 
— O.968U6856 


i|«f 

sinQ 
cos ■* '■ 

Y/R sin/-\ 

X/r cos ^ ^ ^ 

Y/ilsin Q + X/R cos i c ® 


Oo 86601592 
-0.00963126 


0,85168110 

-0.09U5050U 

-O.U38U9138 

-0.89873535 

-0.373U5U82 

O.08U93502 

-0.28851980 


0.81+218232 

-0.1371+1760 






-0.3797U052 
0.00865595 
cos u -0.371081i57 


-0.36928969 

0.12350205 

-O.2I+57876U 


u 


2U8°13’02«,8U5 


253013 «5o« .312 255 ®Uo« 21 ".U 2 U 


u^-U2 = CX 1 
Cfl (radios) 

U3-U1 : O2 

(J'2 (rad^s) 
U2-U1 = O3 
03 (radians) 


2026 8 31" .113 
O.OU262050 


7027 «i 8".579 

0.13011710 


5°008U7"»1+66 

0,0871+9658 



TABI£ 8 
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DETERMINATKaj OF AUXILIARY PARAJaETER, p. 



(T3-T1) 



6 (Ji 03 



Ri ^ cra(2 O ' . O ',) r ,2 cr,(2 cr ,, cr ,) 
60-3 6 cr^ 



Rx 2 

rP 

R^2 



6 f(f. (fsl ^ ^ 

R 22 V2Y(6 CT^CTj) 

^lo (^ 2(2 Cf3- Cfi)/6 0^3 
D2 

(i^l- CT3) 

02(2 01" (/o)/6 (fi 
k yf- (T3^Ti)s *+*+* = 
K (T3-T1) 

«+*+«-/^(T3 -Ti) 3 Yp = 
p 



1.U8058990 

1.50835250 

I.516UU300 

0.0022029it 
O0OO33228I 
0.02237U90 
0oHi8506lii » 

0.19265006 

0,13237266 

0.52U979li8 

O.OU857637 ^ 
0.19731517 
- 0.00225558 
0,25572300 

-0.0017li0U0 «• 

0.195314211 

0.211i3Wi36 

O.9II3U7IO 

0.8305535U 



TABLE 9 
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CCMPUTATICN OF THE ORBITAL PARAMETERS a, w. 



e sin (ui=w) = R-^ (p-Ri)cos 

Rl R3 sin 0 2 



R3 

(P-El) 

cos 




(P-R 3 ) 

R3(P”Ri) ^ 
R3(P-Rx)cos O 2 

Ri(p-R3) 

R3(p-Ri)cos 0 2”^i(P“^3) 

sin Q 2 

R1R3 sin cTp 
e sin (ui - w) 



Io 2311 t 39 li 0 

-O.3862U139 

0.99151*671 

1.216791*93 

- 0 . 1*0088586 

-0.1*7563287 

-0.1i7161221 

-0.1*8779588 

0.01618367 

0.12975025 

0.19l*la897 

0.08321*121 



e cos (ui - w) = (p-Ri)/r1 ~ -0.3171*2521 



c2sin2(ui-^) 


0.00692910 


e^cos^fui-w) 


0.10075876 


e2 rsin2(ui-w)fcos2(ui-iff) 


0 .I 076 B 786 


e >- J 


0,32815828 


* a 


0,93078812 


e sin (ui»^) 


0.08321*121 


sin (ui-tt) 


0.25366177 


e cos (ui-w) 


-0.3171*2520 


cos (ui-w) 


-0,96729298 


(ui-w) 


165018 ' 20 ". 508 


ni 


21 * 8013 '02 ". 81*5 


w 


82051*' 1 * 2 " .337 
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I 
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TABIE 10 



V * u-w 
V/2 



COMPUTATION OF TIME OF EPOCH, To 
Observation #1 Observation #2 

165°18»20"«5o 8 170°19«07".97li 

83°39»10".251i 85°09»33".987 

tan E/2 = ~^l-e)/{l+ e) tan V/2 



Observation #3 

172°li5'39”.087 
86022 •Ii9".51»3 



1-e 
1 + e 

7(l^)/(l+e) 
tan V/2 
tan E/2 
E/2 
E 









7.75557^3 

5.5l59787li 

79°U3*27”.723 

l59°26«55".lUi6 



0.67l81il72 
1.32815828 
0.5058Ui62 
0.7112275U 
11.808i;28 
8 o3981i 7920 
83°12«35".307 
166025* 10". 6lU 






I5»80837li 
U.2li335095 
8U°55*02".65U 
169050*05 " .308 



M = E - e sin E 



E (radians) 
sin E 
e sin E 
M (radians) 
M 



2.782905U 

0.3510I;533 

0.1l5l98ii3 

2.66770698 

153°25*16".711 



2.90ii57022 

0.23li80935 

0.077051i63 

2.82751559 
162000*16" .966 



n - k/V^ 



2.96U7658 
O.I76U868O 
0.05791560 
2.90626098 
166030* 59" .358 



k 




n 



E-e sin E/n-rfc.k.To 0.02237350 
tj (days) 

j = (1,2,3) Dec 18.93692292 



107.0731376 

0.8061^0366 

0.89799981 

119.2351U50 rad. /day 

0.02371378 0 . 02 U 37 li 20 

Dec I8.9382li560 Dec I8.93892li768 



To 



Dec 18.91U5U91 i 2 Dec 18.911453182 Dec I8.91l455056 



TABLE 11 U2 

orbital elements determined frcm published data 



(t-To) - -U.085U6818 days w ® 110° ii8'I;0''.896 

sia 0 . 93 U 7^^25 

(t-To) 2 * 16.69105025 cos w* -0,35529231 

(t-To)3- -68.19075U69 



w = (lU0O.6) + (7°.275){“U.085ii68l8)+(o.000U27)(l6.69105025) 

«(°.22U) (-.35529231) 

w = 110° 57‘ 53”o751i 

Q a (l86o.278)-(U°.896)(-U.0851*68l8)-(°.000298)(l6.69105025) 

-(°.00887) (-.35529281) 

Q . s 206° 16» U3".067 

i = (33°. 210 )-(°,ooo55U)(-U.o85U68l8) + (°,ooooo5) (16,69105025) 

+ (o.00ii73)(.93U75525) 

i = 33° 13’ 00".U02 

e = (.09971)-(.000202)(-li,085U68l8 )+ (o00000003)(16.69105025) 

-.(.oooU5519)(.93U75525) 



e = 0.10011028 

a = q/l-e = 1.053706/.89988972 



U3 

Comparing the above results it Is apparent that the position of the 
orbital plane is correct but that there is a discrepancy in the orienta- 
tion of the orbit in its plane and in the size and shape of the ellipti- 
cal orbit o Reviewing the computations, it can be seen that the auxiliary 
parameter p, which has a value of approximately 0.830, is too small. In 
the case of an elliptical orbit it should be slightly greater than one. 
Since p is a function of the auxiliary parameters uj, it is believed that 
the source of the discrepancy rests here. This would affect the results 
in the computations for and for the argument of the perigee, w, 
also. To avoid continuing the problem with this discrepancy present I 
decided to use the published orbital data from the 82 observations rather 
than to continue with the non-re alistic values for e and a derived from 
the computations. The orbital parameters computed for the times of 
observations in the next chapter are based on those published in Smith- 
sonian Astrophysical Observatory Special Report #78, 
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DETERMINATION OF “fflE GEOGRAPHIC COCKDINATES 
OF THE UNKMCSTO STATION. 



There are several methods available to determine the position of an 
vmkncwn station from observations of a celestial body. Two of these 
methods were devised by William Markowitz of the U. S. Naval Observatory 
for his well known Moon Camera Method and can be adopted for use with 
artificial earth satellites. These methods are called the linear, 
which gives the coordinates of an lonknown station, and the differential, 
which provides a way to find the corrections to an assumed approximate 
position. The differential method has been selected for use in this 
particular problem. As mentioned earlier in this paper, only two 
approximations in the differential method will be made. 

The procedure to determine the coordinates of the vinknovm station 
is as follows: 

1. The orbital elements of the satellite are known for a particular 
time. To, the time of epoch. Now assvuning these elements are correct. 



the orbital elements are then determined for the times of the observa- 
tions from the unknown station. 

2. Using these elements the geocentric right ascension, geocen- 
tric declination, and the range of the satellite from the earth's center, 
will be predicted for the times of the observations from the unknown 
station. 

3« The differences in the computed coordinates and the observed 
coordinates of the satellite will furnish us with a set of coefficients 
for three differential equations to be solved for each observation. The 



solution of the differential equations will produce the increments 




which are added to the approximate coordinates. Using these 










• t > • ^ 



new coordinates, the differential process will be repeated and the resul- 
tant increments will be added to the approximate coordinates again. This 
will be the final result. 

Refer to Figure 6, 

The observer, 0, views the satellite aJ.ong the line OS* at a distance 
of r* The coordinates observed from this position are the apparent right 
ascensicn, (X , and 1iie apparent declination, (5 • PVom the center of 

the earth, E, the satellite would be seen along the line ES* at a range 
of R. The coordinates observed from this position are the geocentric 
right ascension, (X , and the geocentric declination, (5 . The coordi- 
nates of the satellite derived frcsn the orbital elements are the geo- 
centric coordinates and will differ from the observed coordinates determined 
from the camera film by the araotint of arc, SS*, in the orbital plane. 

This difference enables us to determine the observer's position on the 
surface of the earth. This difference between the computed and the 
observed positions of the satellite should be very small but not so 
small that the two positions would be at the zenith of the observer. At 
this position both the computed and observed positions would have the 
same value and the observer's position could not be determined. Another 
indeterminate situation would exist if the satellite was observed at 
the horizon. Here the parallactic angle would be too large and the 
effect of refraction would make the observations useless. [l2, p.3li] » 

Proceeding with the problem, the following data are knowns 

Location of Observer: ^ Curacao, N.W.I. 

Approximate Position: ■ 12° 13' 23 ".21 N 

X » 68° 55'07".86 W 

xa - 0.3519928U2 

ya s 0.91310^082 

Za » 0.210569117 

Explorer I (1958 Alpha) 



Satellite Observed; 
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The approximate rectangular space coordinates, Xg^, are found in 

the same manner as that used to find the coordinates of the known station 
This procedure is discussed in Chapter U and will not be shovm here# 

Using the above data, the orbital parameters for the times of obser- 
vations are ccsnputed. These computations are contained in Table 12» 

With the orbital parameters known, the next step is to compute 
the eccentric and true anomalies of the satellite using this formula 
(U,p.l6lj for the eccentric anomaly: 



where M is in radians » The additional terms in this formula are not 
listed or used in the determination of E« Moulton states that because 
of the anall eccentricity of the orbit the series converges very 
rapidly and the first three terms give an accurate value of Eo 
Then the true anomaly V is computed from the formula: 



The computations of the anomalies for the three observations are listed 
in Table 13 « 

The next step is to determine the geocentric right ascension, geo- 
centric declination, and geocentric distance of the satellite ibr each 
observation using ■Uie formulae: 



E ~ M+e sin 



e.2 



sin 2 M + 



e e o 




siii 0 - sin i sin (w4v) 

a - + arc tan 

R = a (1 - e cos E) 




where a is the semi-major axis of the elliptical orbit in terms of the 
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TABLE 12 



kS 

OBSERVATION 3 



OBSERVATION 1 OBSERVATION 2 



Location 

Time of Observation 
Date of Obsenration 
Smithsonian Date of 
Smithsonian Date of 
(t-To) 
w (Mean) 
sin IT 
cos w 

(t-To)2 

(t-Ta)3 



(U.T.) 

Observation(t) 
Epoch (Tq) 



Curacao,N»W,i. 

23%1®U7®.0U 
10 Dacanbo:* 1$60 
37286.991i29Uia0i 
37291.0 

-ii.005705556 

111°23'27“.72?5 

.93111290 

-.36I473IIO 

16.01*567700 

-6U.27U25751 



!uracao,N.W,l, Curacao ,N.W, I* 

3^51“*5U®.98 23"52“12®.1*50 

18 Decanher I960 
37286.991*58^1*2 

37291.0 

-1*.0()5W58" 

111023' 35" Jt2 00 

.93109929 

-.361*76583 

16.01*332095 
-61*. 26010156 



C 

23 

18 Ifegaitoer I960 
37286.991*38^1*3 
37291.0 

-l*.0056l3657 



1H023'30''.3299 

.93110865 

-.361*71*195 

16.01*1*91*077 
-61*. 26983387 



w » (li*0°.6 ±2)+7°.275(t-To)+0O.000l*27(t-To)^+0°.22l* cos w 



7®.275(t-To) ^ 

0°.000l*27(t-To)2 
OO.22I* cos w 



’^to 



-29^^.11*150792 

0°. 00685150 

-0°o08l69977 
11*00.6 ± 2 
1110.38361*381 
1110 23 ' 01". 118 



-29°.H*083935 
00.00685119 
-O0.08I70220 
ll*0®.6 ± 2 
111°.381*30961* 
111o23»03".^5 



-29°. 13936836 
0°. 00685050 
-00.08170755 
ll*oo,6 ± 2 
1110.38577U59 
U1023'08",789 



Q= (1860.270 ±9) - (-l*'5;896(t-To)-0®.000298(t-To)2+0°.00887 cos w 



l*o.896(t-To) 
00. 000290 (t-To 




COS w 



)2 



19°.6ll93iii*0 

-00.001*78161 
-0O.003235L6 
1860.278 ±9 
2050.88191763 
205°52'5U"^C8 



190.6111*81*1*6 
-0°. 001*78139 
- 0°.00323526 
186°.278 ±9 
205°.88U*6781 
205°52'53".2B1* 



19<^.6lOl*9l*50 

-0°. 001*78091 

-0O.003235U7 
186°.278 ±9 
205°. 8801*7812 
205°52'1*9".721 



i = (33°.210 +3) -0°.00055U(t-To) + 0O,50xl0“^(t-To)2-0^.00l*73 sin w 



0°.000551i(t-To) 
0O.000005(t-TQ)2 
0°.00l*73 sinw 
i To 

H 

it 



00.00221916 
OO.OOOO 8 O 23 
-00.001*1*01*16 
330.210 ± 3 
330.20789523 

33°12»28"ii2S 



00 .00221911 
00 .00008022 

-0O.00l*i*0hll* 
33°.2lO Z 3 
330.20789519 
33°12'28".1^3 



0°. 002 21900 
0°. 00008022 
-0O.00l*l*0lp.0 
330.210 ±3 
330.20789512 
33°12'28".l£2 



e s (,09971±8)-0,000202(t-To)+0.30xl0-7(t-To)2+0.000l*55l9 sinw 



0.000202(t-To) 

0.00000003(t-To)2 
0 . 0001*5519 sin w 
©To 

et 



0.00080915 

0.0000001*8 

0.0001*2383 

0.09971+8 

0.10091*31*6 



0.00080913 

0.0000001*8 

0.0001*2383 

0.09971+0 

0.10091*31*1* 



0.00080909 

0.0000001*8 

0.0001*2383 

0.09971+ 8 
0.10091*31*0 



(continued on page 1*9) 
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TABLE 12 (Continued) 



k9 



OBSERVATION 1 OBSERVATION 2 
M ■ (.9U0l6) + (13.U35l3ii)(t^To>f0.000283?(t-To)2 -0.11jcl0“^(t-To)' 



(13.U35l3U)(t-To), 
(0.0002839)(t-To)^ 
(0.0000011) (t-To)3 

(0.00000063) cos w 

MTo 

M 

Mt (Revolutions) 

Mt 

Mt 



-53.81719091 

o.ooii55537 

0.00007070 

0,00000023 

0.9U016 

-53.81256U61 

0.12759539 

U5°.93U3U0U0 

U5®56'03'*jS25 



“53.81595623 
o.ooii555l6 
0,00007070 
0,00000023 
0.9UOI6 
-53.811330m 
0.12882986- 
U6°.3787ii960 
1^6022 »U3''J69 



observation 3 

’-0.000631x10“^ 
cos w 

-53.81323966 

o,ooi;551;70 

0.00007069 

0,00000023 

0.9U016 

-53.8o86lUoii 

0.1315U596 

It7°.35651i560 

U7°21»23".55ii 
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TABLE 13 50 

COMPUTATICN OF ECCENTRIC AND TRUE ANOMALIES 



E = 


M + e 


OBSERVATION 1 
^e2 

sin sin 2 M 

0.01018958 


OBSERVATION 2 
0.01018958 


OBSERVATION 3 
0.01018957 


e2/2 




0.00509U79 


0.005091*79 


0 . 005091*79 


1 + e 




I.l009ii3lt6 


1.10091*314* 


1 , 10091 * 31*0 


1 - e 




0.8990565U 


0.89905656 


0.89905660 


(l+e)/(l-e) , 

(l+e)/(l-e)]2 
'■ sin M ■' 




1.22^5510-9 


1.221*551*11* 


1 . 221 * 551 * 01 * 




1.10659577 


1.10659571* 


1.10659570 




0.7185U327 


0 . 72391601 * 


0.73558351 


sin 2M 




0.999U6819 


0.99881*211 


0.99661861* 


e sin M 




O. 0725322 U 


0 . 073071*58 


0 . 071*25230 


(e2/2)sln 2 M 




0.00509208 


0.00508889 


0.00507756 


E (radians) 




0.87932980 


0.88762535 


0.90585750 


E 




50938I8863U 


50985718631* 


51990181159 


£ 




5Q®22«5U".79L 


50O5l'25''.8?l 


51°51*>o6'».522 


E/2 




25°U'27".396 


25°25'1*2''.$66 


25®57’03".261 


^ tan ^/Z a 1 

[(!■+« )/(l-e)l 2tanE/2*tan 

L v/2 J 




O 0 U 7037122 


0 . 1 * 7514*692 


0 . 1 * 8667231 * 


V/2 


0.52051080 


0 . 52612751 * 


0.53851*952 




270291 5o'»,872 


27°1*5'00".338 


28018' 16«.789 


V 




51i059'Ul".7Wi 


55°30«00".676 


56o36'33”.578 



■»K- 
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CQMPUTATICN OF GEOCENTRIC DECUNATICN, GEOCENTRIC RIGHT 
ASCENSIOJ, iND GEOCENTRIC DISTANCE OF THE SATELLITE. 



OBSERVATTCN 1 OBSERVATION 2 OBSERVATION 3 
sin (5 ~ sin i sin (w+V ) 



(w+V) 


166°22'1|2".862 


I66°53'0l|".191 


167059 'U 2" .367 


sin (w+v ) 


0.23550558 


O. 22691 U 83 


0.20799531 


sin i 


0.51i767853 


0.51i767853 


0.5Ii767853 




0.12898135 


O.I 2 U 27638 


0.1139Hi57 


6 


702U'38".636 


7008 ' 20". 290 


6°32»27".682 


a = 


Q+ arctan cos i tan (w+v )] 




tan (w+v ) 


-O.2I1232136 


-0.2329925U 


-0.21261i591 


cos i 


0.83668886 


0.83668886 


0.83668886 


cos i tan (w+v ) 


-0.2027U758 


-0.19U91|226 


-0.1779181|6 


Q r 1 


1 205°52>51|".903 


205052' 53 ". 281 | 


205°52'U9".721 


arctan cos i tan(w+v ) 


168032' 19". 601 


168®58'08".322 


169 ® 51i'lA". 720 


a '■ 


liiO25'li|H,50l| 
a = q/ (1-e) 


lUo 5 l' 01".606 


l5°li7'31''J4lil 


q (Ik) 


6.72127351 


6 . 721273 U 9 


6.721273U9 


q (earth's radius) 


1.05375739 


1.05375739 


1.05375738 


a 


1.17207021 


1.17207018 


I0I7207OI2 




R ■ a (1 -» e cos E) 






e 


0.1009li3l;6 


0.1009U3liU 


0.1009U3UO 


E 


50®22'5Ii".791 


500511 25".871 


51°5U'06".522 


cos £ 


O. 6376675 U 


0.63125552 


0.61701099 


e cos E 


0.061i36837 


O. 06372 UO 


0.06228319 


1 - e cos E 


0.93563163 


0.93627890 


0.93771681 


R 


1.09662596 


1.097381i58 


1.09906985 



52 



earth's radius. The value «a" can be computed from the formula: 

a = q/(l-e) 

where q is "the distance of the satellite at perigee. The computation 
of these values are shovm in Table lit. the Greenwich hour angle 

of the vernal equinox j is tabulated in Table 15. The computations of 
the sidereal times of observations from which this value is derived is 
not shown or listed. The sidereal time was established in the same 
manner as it wais for the observations from the knovm station. 

Now the counted apparent right ascension and declination must 

be found using the formulae [l2, p.B?] : 

4 .^ fry' ry ^ ® Xa sin G cos G 

an ( OC c- Ca ) _x^ cos (i+ya sin G+fe cos ^ 

and ^ 

tan 6c = - 0^ ) m 6 

j ' -ix^ cos G + sin G+R cos ^ 

Then the differeoces between the con^)uted and the observed values 

are determined byr 



A Of = a - a. 
A 6' = 6'- 6c 



The values and computations of 0 q and OC ^ and the increments 
A CX and ^ 6 are tabulated in Tables 15 and l6. 

The differential equations to be solved [l25P«36j are: 

I 

ai dXa -t a2 dja = a3 dQC ^ 

b]^ dx^ + b2 dya+^3 = b^ dC^ + b^ d6 
The coefficients of these equations are; 

a]^ ® sin G + tan ( CX - (X ) cos G 

I 

a 2 = cos G - tan ( OC “ 00 sin G 

a 3 = sec2 (OC-OC )(R cos (5 - cos G + sin G) 
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- tan (5 cos G 
b2 - - tan (5 sin 0- 
b3 = - cos ( CX - CX ) 

- sin (CX-CX) (R sin 6 - z^) 

b^ = sec^ (5 (R cos S -x^ cos G + sin G) 

The computation of these coefficients are contained in Table 16. The 
solution of the first approximation differential equations is shovm in 
Table 17. AGIO Technical Report No. 86 suggests these equatirais be 
solved by Doolittle’s method but being more familiar with the Gaussian 
solution of nomal equations, I decided to use the latter method. The 
solution of the normal and reduced equations yields the increments 
Axaa Aya» snd Az^. These increments are then added to the origin- 
al approximate coordinates to find the corrected position of the obseiver. 
The corrected space coordinates areg 

Xo yo 2o 

+ 0.351U62753 + 0.912606236 + 0.21092082U 

Since the International Ellipsoid was used as the reference ellip- 
soid, the xa, yaj Za coordinates can be multiplied by 6,378,388 meters, 
the semi-major axis of the ellipsoid, to find the metric distaice of P, 
the radius of the observer's position fl*om the earth's center. This has 
not been done because it is of little interest in this problem. How- 
ever, the distance of the observer fi*om the earth's center is found in 
units of the earth's radius by the equations 

P ~~\Jxl -f y^ + zg' » 1.000U31808 

Then the jresultant geographic coordinates of the observer may be easily 
derived as follows: 

sin (j)" = ^°/p = + 0.21082979 
tan 0' = +0.21^67762 



* 



f 






TABIE 

/(APPAEENT RIGHT ASCENSION, OC cl 
SOLUTICN FCR COMPUTED dECLINATICN, (5|. 



FIRST APPROXIMATICN 






CC 

6 

R , 



9g (GHA Y ) 



OBSERVATION 1 

7°2U«38».636 
1.09662^96 
3^50 U6«30". 
-30° ^I'OO". 
8^°31'5U”.97? 



OBSERVATION 2 

lU°5l'01".6o6 
7008' 20" .290 

1.09738Ji^8 
3^9° 39' 00". 
-32° ^2' 00". 
85033 '51i".l|05 



OBSERVATION 3 

i5oii7'3r'.liU 
6032* 27". 682 
1.09906985 
8° 03' 00". 
-36° 39' 00". 

85°38'17'M6o 



cx- g= G 

sin G 


-71006' U0«.U71 


-7001;2'52".799 


-69050 'ii5“ *719 


-0.9U6lIi889 


-0.91;388553 


-0.93877031; 


cos G 


0.32373178/ 


0.33027279 


0.3Ui;5i;l;07 


slnr 
COS ^ 

( Gc'- a ) 


0.12898135 


O.I2U27638 


0.113911;57 


0.99161^703 


0.99221;765 


0.993U9055 


3iil021' 15“ .U96 


3Ui;01i7»58".39lt 


-7°Ul;'3rMUa 


sin (a' -a) 


-0.31971509 


-0.26219670 


-0.13i;7l37i; 


cos (C(' -00 


O.9U751372 


0.965011;95 


0.99088U56 


tan (a; -a) 


-0.337li2529 


-0.27170236 


-0.13595301 


sec2(a' -a) 


1.11385581; 


1.07382217 


1.0181;8322 


tan O' 


-0.59730303 


-0.61;6l0U06 


-0.7Ul;020l;0 


s®c^ 


1.35677091 


1.U17U501;6 


1.55356636 



tan (a« -a) 



sin G-fy^ cos G 
=Xa cos 0 + 7a sin G+R cos ^ 





>. 


+ O.3519928U2 
+ 0,913105082 
+ 0.210569117 


^ 






^ 


^ 




^ 




x^ sin G 
7 a cos G 
-Xa cos G 
7 a sin G 

R cos Q 


- 0 . 33303761 ; 

0.29560113 

-0.11395127 

-0.86393336 

3 .. 087 U 6588 


- 0 . 3322 i ;095 

0.30157376 

-0.11625366 

-0.86186667 

1.08887727 


- 0 . 330 l;l ;037 

0 . 311 ; 60 l; 9 l; 

-0.12127705 

- 0.85719579 

1.09191551 


XasinG+7a cos G 
-XaC osGf7asinG+ 
R cos (5 


- 0 . 037 i; 365 l 

0.10958125 


-0.03066719 
O.IIO7569I; " 


- 0 . 015835 U 3 

0 . 113 U ;267 


tm ( dc - OC ) 

i(x' - a) 

«'c 


- 0 . 31 A 6321 U; 

3 la° 08 « 17 ".l ;06 

355 ° 33 ’ 31”.910 


-0.27688730 

314 ;° 31 ' 23".757 

359022! 25" .363 


- 0.13958972 
-7056' 47" .593 
7050' 43". 848 




tan 6' = 00s 


(CX'c - OC ) (R sin 


6 - za) 



-Xa cos G+7a sin G+R cos 



TABLE 1$ (Continued) 
FIRST APPROXIMATION 
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OBSERVATION 1 OBSERVATION 2 



cos (Cl'c-CX) 0.9U630093 

R sin Q 0«lltl)|l|lt3Q 

(R sin (5 -Za) -0.06912U82 

-Xg, COS G + yg sin G + 

R cos 5 0.10958125 

cos( OC' fOC)^ 

(R sin O-Zg) -0.065U1288 



0796373889 

0.13637898 

-O. 07 U 9014 

0 . 11075691 * 

-0.0711*9992 



6'c 




- 0 . 596931*97 -0.61*555702 

- 30 O 50 ' 01*" .037 -32°50' 1*0" .376 



OBSERVATION 3 

0.9903971*6 

0.12520007 

-0.08536905 

0.1131*1*267 

-0.081*51*929 



-0.71*530U13 

-36°1*1'50".335 



TABIE 16 
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CCMPUTATIQN OF COEFFICIENTS OF DIFFERENTIAL EQUATIONS 
FIRST APPROXIMATION 





OBSERVATION 1 


OBSERVATION 2 


OBSERVATION 3 


^^^(radians) 


ax ■ sin G+tan iOC ' - OC ) cos G 




0 ° 13 ' 58".090 
0°. 23280278 
O.OOU06318 
- 0 ° 01 * 0 U '«.037 
-00.01778806 

- 0 , 000310 li 6 


0 ° 16 ' 3 U ".637 

0°. 27628805 

0 . 00 U 8221 U 
- 0 ° 01 « 19 ”. 62 U 
-0°. 02211778 
-0.00038603 


0 ° 12 * 16",152 

00 . 20 Ui 8666 

0.00356897 
0002 ' 50 '' .335 
O0.0U731527 
0.00082581 


sin G . 

tan(OC-OOcos G 


-O.9U6H4889 

-0.10923529 

- 1 . 05538 U 18 


- 0 . 9 U 388553 

-0.08973590 

-I.O3362H43 


- 0 . 93877011 t 

-0.0U68U180 

- 0 . 9856119 U 




32 = cos G - 


tan ( (X - OC ) sin G 




cos G , 

tan(CX-OOsin G 
a .2 


0.32373178 

0 . 31925 U 56 

O.OOUi 7722 


0.33027279 

0 . 256 U 5593 

0.07381686 


0 . 3 UU 5 UU 07 

0.12762863 

0 . 216915 UU 




33 = sec 2 (C( 


*-C 0 (R cos( 5 “Xa cos G+7a sin G 


sec 2 (n( -CO 

R cos^ 

Xg ^ COS G 
7 a sin G 

Rcosn-Xa COS G + 
7 a sin G 


1 . 1138558 U 

1 . 087 U 6588 

0.11395127 

-0.86393336 

0.10958125 

0.12205772 


1.07382217 

1.08887727 

O.H625366 

-0.86186667 

0 . 1107569 U 

0.11893326 


1 . 0 l 81 i 8322 

1.09191551 

0.12127705 

- 0.85717579 

0 . 113 Uii 267 

0 ,U 55391 a 6 




bx .= tan (5 ' 


cos G 






-0,19336597 


-0.21339059 


- 0 . 2563 U 782 




b2 = - tan (5 


t 

sin G 




bg 


-0.56513760 


-O.6098U827 


- 0 . 698 U 6 Um 




b3 = - cos ( a ‘-( X ) 






- 0 . 9 li 75 l 372 


- 0 . 9650 lUli 5 


- 0 . 99088 U 56 




bl^ = sin (CX 


-CX )(R sin 6 -Za) 




sin (CX'-OO 
R sin Q 

Za r 

R sin 0 -Za 
bU 


-0.31971509 

0 ,llillM 30 

0.21056912 

-O.06912I482 

0.02210025 


-0.26219670 

0.13637898 

0.21056912 

- 0 . 07 Ul 901 ii 

0 . 019 U 52 U 1 


- 0 . 13 U 7137 U 

0,12520007 

0.21056912 

-0,08536905 

0.01150038 



1 



« 
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TABIE 16 (Continued) 



FIRST APFROXIMATICN 



OBSERVATION 1 OBSERVATION 2 

= sec^ (5^ (R cos (5 - cos G+y^ sin 

1.35677091 l.la7U5oU6 

R cos (5 - x„ cos G+ 

7a sin G 0.10958125 0.11075691 

O.IU867665 O.156992U8 



OBSERVATION 3 
G) 

1.55356636 

0.113iiU267 

0.17621072 



f 



I 



{ 



TABIE 17 
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SOLUTION OF DIFFERENTIAL EQUATIONS 
FIRST APPROXIMATICN 



b^Aa' 

'=UAa'+*’5A(5' 



OBSERVATICN 1 

O.OOOU 959 ii 

0.00008980 

-O.OOOOU 616 

O.OOOOU 36 U 



OBSERVATION 2 

0.00057351 

0,00009380 

-0.00006060 

0.00003320 



OBSERVATION 3 
O.OOOiil 236 

o.ooooUloU 

0 . 0001 U 55 U 

0.000186^8 



OBSERVATICN EQUATIONS 

aidxa + ag^ya = ^ 

bidxa + b2dya + b3 dza == b|^ d CX' 



b^d 6 ^ 



dxq 



dxa 


8 ya 


dza 


1 


s 


-I.05538UI8 

- 0.19336597 

-I.03362IU3 

-0.21339059 

- 0 . 9856119 U 

-O.2563U782 


O.OOlUi 7722 

- 0.56513760 

0.07381686 

- 0 . 6098 ii 827 

0 , 2 l 69 l 5 Ui 

- 0 . 698 U 61 ilii 


- 0 . 9 ii 75 l 372 

- 0 . 965011 i 1;5 

- 0 . 99088 U 56 


-O.OOOU 959 U 

-0.00001436U 

-0.00057351 

-0.00003320 

-O.OOOiil 236 

- 0,00018658 


-1, 051^0290 
-1.70606093 
-0.96037808 
-1.78828651 
-0.76910886 
- 1 . 9 U 588310 






NORMAL EQUATIONS" 




dXa 


dya 


dZa 


z 


E 


3.302280071 


0 . 1236 U 5699 

1.231668857 


O.6U3153009 

1.816085355 

2.810887350 


0.001585979 

O.OOOOU1227 

0.000258267 

0.000000783 


it. 07066 U 758 
3 . 171140137 
5.270383981 
0.001886255 






REDUCED EQUATIONS 






dya 


dZa 




E 


3.302280071 

- 0 ;i 9 it 760285 

-0.000U80268 


0 . 1236 U 5699 

1.227039250 

-I.U 6 oU 2930 ^ 

a.ooooiU 797 


0 . 6 U 3153009 

I.79200U08U 

0 . 068531 ii 01 

0.000351707 


0.001585979 

-0.000018156 

- 0 . 00002 lil 03 

0.000000013 


U. 07066 U 758 

3.019025178 

0.068507298 

0.000000012 


dxa 

-0.000530089 


dya 

- 0 , 000 U 9881 i 6 


dZa 

+ 0.000351707 







Xa O.3519928U2 
Axa -0.000530089 

xo + O.351U62753 



Ayt 



0.913105082 

-O.OOOU 9881 i 6 
7o+ 0.912606236 



a; 



'a 



0.210569117 

0,000351707 



Zo+ 0.21092082li 



( 



« 

« 



r 

k 
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tan tan <^/l-e2 = 0«2l567762/0.99327733 = 0.21713736 

tan A = 7o/xo = 2o59659l4U6 

which result ins 

<j) - 12° 15' 03'M68 N 

A = 68° 56' m«o336 W 

The differential method was repeated again for a second approxijna- 
tion. The corrected space coordinates from the first differential 
process were used as the approximate obseirver's position. The data 
used in and the computations of the second differential approximation 
are listed in Tables 18 ^ 19, and 20. 

The corrected and final results of the second approximation are: 

Xo Jq Zo 

+ 0.35lli8350U + 0,912709058 + 0.2108U9193 

P = + yo^ + zo^' = 1.000517795 

sin = ^o/p - 0.2107U007 

tan 0'= 0.21558157 

tan 0 = tan 0/l-e2 = 0.21558157/0.99327733 = 0.2170U066 
tanA = Wxo = 2.59673370 

The corrected and final geographic coordinates are; 

0 = 12° lU' Ui".H8 N 
A = 68° 56« 18".0U5 W 

Further approximations would probably yield changes to these 
coordinates but as stated previously the assumption that two approxi- 
mations are sufficient has been made. 

As a note of interest, and to satisfy my own curiosity, further 
calculations were made using the linear method but are not included in 
this paper. First, I used ngr own computed value of the Gaussian gravi- 
tational constant, k, in the formula to determine "a": 



I 



I 



i 

f 



I 



1 
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TABLE 18 



o^TTTm-r^T nr. /^APPARENT RIGHT ASCENSION, CXc 

SOLUTION OF COMPUTED DECLINATION, (5'^ 



SECOND APPROXIMATION 



OBSERVATION 1 OBSERVATIOI 2 OBSERVATION 3 



tan ( C(c-C( ) _ xa sin G -t- ya cos ,G 







“Xa cos G + 7a sin G + 


R cos^ 









+ 0 . 3511*62753 




7a 






+ 0.912606236 


<: 


za 




^ 


+ 0 . 210920821 * 





Xj^sln G 




“0.33253609 


- 0 . 33171*060 


- 0 . 32991 * 271 * 


7aC0s G 




0.295U3961t 


0 . 3011*0901 


0 . 311 * 1*3307 


-Xa cos G 




-0.11377966 


-0.11607858 


- 0 . 121091 * 1*1 


7a sin G 




-O. 863 I 16138 


-0.86139582 


-0.8567271*9 


R cos ^ 




1,0871*6588 


1,08887727 


1.09191551 


x^sinG + 7 aC 0 sG 




-0.0370961*5 


- 0.03033159 


- 0,01550967 


-XaCOsG-^a sin 


G + 








R cos Q 


0.110221*81* 


0 . 1111*0287 


0 , 111*09361 


tan (a'c-OO 




“0.33655251* 


- 0.27226938 


- 0,13593811 


(a>,-oo 


f 


3ia°23>57".123 


3l*l*°l*6»09”.l*92 


-7°1*1*»28" Jb2i* 




«c 


355°1*9«11",627 


359°37»11''J098 


80 03«03"JD17 




tan (5^ 


_ cos{Q(!^(X )(R sin 6 








-Xa cos G + 7a 


sin G+R cos 




cos (OC'c-OO 




0 . 91*776396 


0.961*87589 


0.99088653 


R sin Q 




0.11*11*1*1*30 


0.13637898 


0.12520007 


R sinS ~Za 




- 0 . 0691*7652 


-0.071*51*181* 


-0.08572075 


cos((Xl“Ou(R sinO-Za) 


-0.06581*731* 


- 0.07192362 


-0.081*93951* 


tan 


rt 


- 0.59739111 


-0.61*561730 


- 0 . 71 * 1 * 1*7237 


L/ C 


6c 


-30O51«13".391 


-32°50» 1*9" JL52 


-36°39«59".995 



Pi 



TABLE 19 
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CCMPQTATION OF COEFFICIENTS OF DIFFERENTIAL EQUATIONS 
SECOND APPROXIMATION 





OBSERVATICN 1 
a^^ = sin G + tan ( CX'- 


OBSERVATICN 2 
"(X) cos G 


OBSERVATION 3 


a-a'c=4a' 


-0°02»U1"o627 


0°01«1+8''.912 


-0°00'03".017 


« 1 ACC' 

Act (radians) 


-OO.OI+I+8I+639 


0°, 03025333 


-0® .00083806 


-0.00078272 


0.00052802 


-0.000011+63 


6'-&-A<5; 
,, A<5' 


0<^00«13''.391 


-0O01«10'».8l+8 


0O00»59".995 


0°.00371972 


-0°. 01968000 


0°. 01666528 


(radians) 


0. 000061+92 


-0.00031+381 


0.00029086 


al 


-1.055381+18 

a2 = cos G - tan (CX 


-1.0336211+3 
^-CX) sin G 


-0.98561191+ 


®2 


O.Oi+1+1+7722 


0.07381686 


0,2169151+1+ 



a3 = sec2 (CX'-‘(X){R cos 6 - cos G + ya sin G 



sec ^ ( a -CO 
R cos (5 
Xa cos G 
y sin G 

(k cos (5“Xa cos G+ 
ya sin G) 



1.1138^58ii 1.07382217 

1.087U6588 1.08887727 

0.11377966 0.11607858 

=Oo863U6138 =0.86139582 

0.11022U8U 0.111it0287 

0.12277U58 0.11962687 

= tan S cos G 



1 . 0181+8322 

1.09191551 

O.I2IO9U+I 

-0.8567271+9 

0 . 111+09361 

0.1162021+3 



b 



1 



-0.19336597 -0.21339059 -0.25631+782 

b2 =-tan S sin G 



b2 



-0.56513760 =0.60981+827 -O.698I+6I+H+ 

b^ = - cos (O^^-OO 



^3 



-0.91+751372 -0.965011+1+5 -0.990881+56 

b|^ = sin (CX^”00 (R sin 6 



sin ( a'-oo 

R sinQ 

sin (5 - Za 
bl+ 



- 0.31971509 

O.ll+li+i+i+30 

0.21092082 

- 0 . 0691+7652 

0.02221269 



-0.26219670 

0.13637898 

0.21092082 

-O.07I+5I+I8I+ 

0.01951+1+62 



- 0 . 131 + 71371 + 

0,12520007 

0.21092082 

-0.08572075 

0.01151+776 



TABLE 19 (Continued) 
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OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 

= sec^ (5\r cos (5~3Ca cos G+y^ sin G) 
1 o 35677091 l*ia7U50li6 lo55356636 

(R cos 0 "Xa cos G+ 

73 sin G) 0,11022U8U 0.111U0287 0,111;09361 

O.Ht95U986 0.1^790805 0.17725199 



TABIE 20 
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SOLUTICN OF DIFFERENTIAL EQUATIONS 



SECOND APPROXIMATION 



OBSERVATION 1 OBSERVATION 



OBSERVATION 3 



. I 



-0.00009610 

-0.00001739 

0.00000971 

-0.00000768 



0.00006317 

0.00001032 

- 0 . 000051 i 29 

-O.OOOOi ^397 



OBSERVATICN EQUATIONS 



-0.00000170 

-0.00000170 

0.00005155 

0,00005138 



al + a2 dya = a3 dQ:' . 

bi dXa + b2 dy^^ + dZa = b|^ dCX + b^ dO 

dxa dya dz^ i S 



-I.05538UI8 

- 0.19336597 

- 1 . 033621 ]t 3 

-0.21339059 

- 0 . 9856 ll 9 li 

- 0 . 2563 U 782 



0,00iUi7722 

- 0.56513760 

0.07381686 

- 0 . 6098 ii 827 

0 . 2 l 69 l 51 *ii 

- 0 . 698 U 6 hlli 



- 0 . 9 U 751372 

- 0 . 96501 UU 5 

- 0 , 99088 ii 56 



0.00009610 
0. 00000768 
-0.00006317 
O.OOOOU397 
0.00000170 
-0.00005138 



-1.05081086 

-1,70600961 

- 0 . 9598677 U 

-I.7882093U 

-0.76869U80 

- 1 . 9 U 57 U 790 



dXa 



NORMAL EQUATICNS 
dya dZa 1 



Z 



3.302280071 



0 .l 2361 i 5699 

1.231668857 



0 . 6 U 3153009 

1.816085355 

2.810887350 



-0.000035501 

0.000000868 

0.000001203 

0.000000018 



U.O690U3279 

3 . 1711 t 00778 

5.270126916 

-O.OOOO33UI2 



dx 



a 



REDUCED EQUATIONS 
dya da a 



dXa 3 . 302280071 
dya - 0 . 037 UU 2523 
dza - 0 . 19 U 760285 
0.000010750 



0 . 1236 U 5699 

1.227039250 

-I.it 60 U 2936 $ 

-0.000001790 



0 . 61 i 3 l 53009 

1 . 79200 U 081 i 

0.0685311i01 

- 0 . 0000^631 

dz. 



i 

-0.000035501 

0.000002197 

0.000001^909 

0.000000017 



li. 0690 U 3279 

3 . 0190 U 5531 

0.068536309 

0.000000016 



dXa dya —a 

+ 0.000020851 + 0.000102822 -0.000071631 

Xa 0.3511*62753 ya 0.912606236 Za 0.210920821* 
Axa + 0.000020851 0.000102822 -0.000071631 

xo + 0.3511*83501* yo +0.912709058 Zo + 0.2108U9193 




where n is the velocity of the mean anomaly^ M» The resultant "a" from 
this formula has the value lol7190297 in contrast with the value 
1*17207018 which was computed from the foimila a =q/l-e. The final 
result of geographic coordinates using this value and iiie linear method 
is listed in the stunmary* 

Second, the linear nethod was used again in which the eccentricity 
derived from the mean orbital elements was used instead of the smooth 
value for e« The mean eccentricity for the time of observation is 
0,10067U97« The resultant coordinates from this experiment are also 
listed in the summary* 

Summary: 



Method 




Approximate 
Differential (1st) 
Differential (2nd) 



Linear (using computed k) 12^ 15® 06''*35l 
Linear (with mean e) 12° 15® U5"*2?6 



12° 13® 23”. 21 N 68° 55® 07". 86 

12° 15 » 03 ".168 68° 56® 1U".336 
12° Hi® Ui"*ll8 68° 56® 18»«*0!;5 
12° 15® 06''. 351 68° 56® 11«.075 
12° 15® U5".2?6 68° 57’ 25". 278 
12° 05« 50",Ii 68° 50® 1U".0 



Actual 



6o THE INTER-CCWTIMENTAL TIE 



The last step of this thesis is to effect the “tte" between the 
two stations, the geographic coordinates of which are now known o This 
problem can be considered to be one of geometric geodesy or what is 
more cemnnonly known as the "determination of long lines" and is considered 
to be a part of higher geodesy® 

H. F. Rainsford states there are two main problems in the field of 
long lines [6, p.l2 s 

a) Given the geographical coordinates of two points, 
to find the distance and azimuths between themj 
commonly called the reverse (or inverse) problemj 

b) Given the geographical coordinates of one point 
and the distance and 'azimuth from it to another, 
to find the coordinates of the second pointj 
commonly called the direct problem® 

Since only the geographic coordinates of the two points are known 
in this case, the reverse problem applies® 

The following data are givens 



back azimuth, CC 21 , at Point 2j and determine the distance, s, between 
the two points® 

Since the reference surface chosen for this paper is the Inter- 
national Ellipsoid, the distance, s, between the points will be one of 
the basic geodetic curves called the geodesic® It is defined as [l3,p®2]s 

a line cn any mathematically defined geodetic reference 
surface that defines the shortest distance between any 
two points on that reference surface® It is geometrically 
equivalent to a great circle on the sphere. It is a unique 
line on the reference ellipsoid that ccntains the normal 
to the ellipsoid at each point® 



Point 1 



Point 2 




Q,.,. . Africa 




'^1 

Problem; 



^ ' 

Find the forward azimuth, ^12> at Point Ij find the 
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Refer to Figure 7 » It can be seen that the line, s, connecting the 
two points. Pi and P25 contains the normal to the ellipsoid along its 
entirety* The forwai*d and back azimuths, CX 12 and 21 respectively, 
are shown at the points* The curvature of the geodesic in the diagram 
is exaggerated for clarity of figure only* 




Figure 7 o 

It should be noted that the distance difference between a geodesic 
and a plane curve is very small but the difference in the azimuths of 
the two would be appreciable* As the distance between points increases, 
the difference between the two lines (geodesic and plane cui^re) diverges 
rapidly* 

There are several methods available in geometric geodesy to deter~ 
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mine the geodesic and the azimuths but because this particular problem 
is concerned with the inverse problem and the distance between the two 
points is approximately 6^000 miles^ then a choice of only two methods 
is available. These are the Hainsford Method and Sodano’s Fousrbh 
Ifethod. Another method that solves for distances greater than ^,000 
miles is Rudoe’s Method but since this procedure treats the normal 
section only, it is not applicable for this problem. The method 
chosen was Sodano’s Fourth which is a non-iterative solution. The 
accuracy is of the order of i0”.l in azimuth and i 0.5 meters at 
6,000 miles o 





Figure 8. 

Proceeding with the solution, an atixiliary spherical triangle is 
employed. In our problem, the westward point. Pi, is Curacao and the 
eastward point, P 2 is Olifantsfontein, South Africa. The difference in 
longitude, L, is always positive, ^13, p.Ulj . 

The first step is to find the parametric latitudes, y^and /?2, 



using the formula 



tan y3 = (1 - f) tan (f) 
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with the results 

1 = 12° 07“ 33''.?21 

^2 = -25° 53* 01'M89 

Using the parametric latitvides, the approximate spherical distance, Q , 
is found from [l3,p«lil] % 

cos 0 — sin fi, + cos cos /^2. COS L = -Oo 2017305506 



and 



% 



A A 



= 0.77515W582 



>Ap 



" = 0 ° 18 » 20”,866 



sin Qjz = cosfi} cos sin^ ^ sin ^ 

The resulting value for ^ is 101° 38* 17"«336» 

Now A must be found using the fonnula A^I*+( A “I*) where [l3,p»U2j ; 

( A -1')= ^ sin 0 4- (F-^E) G 

K2 

where P = 2062614.", 80625 » The value for A is: 

L = 97° 11’ 09'MU5 

{\_j~X) - 0° 18 « 20". 866 

X = 970 29s 30 "o 0 H 

Then the forward azimuth OC12 is computed from the formula 

[13, P.U2] 8 

cot CX 12 ~ [tan cos - cos A sin/^/.^sin A = 11U°16®06",607 
The back azimuth is [l3s p,li2 ] t 

cot GC 21 [sin /?2 cos A - cos A tan ] ^sin A = 277°U9’56",5oU 
The next step is to find the true spherical distance, 0©, using A~ 
97°29’30".011 instead of L = 97 °H* 09 ",11:5 in the same formula as 



before : 



cos Oo sxnJ3j son Jj.^ + cos cosj^z cos A = -0,2063866912 
Oo - 101° 51:* 38",i:72 

The final step is to find the geodesic distance, s, from the formula 

[13, p»l:3] 8 
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s =b(Ao Q o'^ Bq sin Qq cos 2 CT - Cq sin 2 cos h Cf ) 
s = 11,309,3U2o6228 meters = 6,102<,ltlt21 nautical miles. 
Refer to the following tables for the detailed computations of the 
constants, coefficients, and the unknorms. 

Summarizing the results: 



Point #1 



Point #2 




acao, N.W.i, 

12° lii'UU".ll8 N 
68° 56'l8”.Oii5 W 



Olifantsfontein, South Africa 
Gi25° 57 • 3U''.700 S 
28° lii' 51". 100 E 



CX 12 = HU° 16' 06".607 
CC 21 = 277° 56".5oit 

s = 11,309, 3U2 .6228 meters. 



This concludes the computations. 



cr p 



TABLE 21 
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LIST OF CONSTANTS 

FOR USE WITH SQDANO»S FOURTH METHCD 



INTERNATIONAL ELLIPSOID 



l/a= 297.0 

= 6,378,388.0000 meters 

= 6,356,911.9U62 meters 

e = 0.0819918898 e2 = 0.0067226700 

e' = 0.0822688896 e' 2 = 0.0067681702 



TT, _ 16 N - ^ 2 

l6eW^ 



237.2388918 





= 70519.511U5 


e2(l6e^2 + e^^ 


K-. _ I6e2 n2+6 2 

‘‘s -5-2 


= ii.9865206it9 


2 e 

- ife5Ti2+e'''2~ 


= O.UOIO8I263O 


jr^ _ l6e2 n2 

l6e%2 + e' ^ 


= 0.799U593686 


f _ 16 .2 h2 


= 3.9865206U9 


6 - 7-5 



e 



TABLE 22 



COMPUTATION OF VARIOUS TERMS 
J^l cos ^2 



A = 



COS/ 



'2. sin L 



son 



cos/^l 
sin L 



cos 



A 



e 



cosy 9 j cos/^r sin L , 

(cos Pi cos Rii sin L)/sin Q 



B = a 2 = 
C = (cos 0 



0,8796086600 

0.9921li55600 

0.8726998266 

O.89IOI82190 

O. 793913 U 666 

B cos 0 

U.9865206U9 

-0.1601^66007 

-O.Ola 5739 U 99 

-0.0083372662 



K3 

B cos^ ^ 

(cos^ - B cos u h 
(cos 0 -B cos 0 K3 

D= -Kl| sin sin^j, 

K, - ^ O.UOIO812630 

sHn a, sin pi -0.0917010172 

s±ci^l siji^2, ~ 0.0367797603 

E = -Kt^ sin A s±D.P% 

0 . 799 U 593686 

0.0917010172 



s^ Pi sin P% 

-K 5 sixiJ0, siny9^ 



G 

K6 

CK6 



(radians) 
0 /sin0 



F = CK6 



G =^Vsin^ 



0.0733116370 



-0.0083372662 

3 . 98602 O 6 U 9 

-0.0332366839 



101038 ' 17 " .336 

I. 7739203 UU 

3 .IU 6793387 

3 . 2128 U 61 O 0 



TABLE 23 
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CQMPUTATIQM 



f 

(i-fX 

tan m 



PARAMETRIC UTITUDE & APPROXIMATE SPHERICAL DISTANCE 
Parametric Latitude ( 
tan ^ = (1-f) tan cf) 

0.00336700311 
0.9966329966 

— O.U8686o87 
_ - 0.^852216078 
^^=- 25 ° 53 '01" .189 



0,21558157 
0.21U8557061 
P = 12° 07 '33 ".521 



L=AA 

sin L 
cos L 
sin 
cos 



APPROXIMATE SPHERICAL DISTANCE 
cos 0 — sin sin U-> + 






0.2100618662 

0.9776880956 



2^ T COS cos /^iCOS 

97° H'09".1U5 
0,99214556 
0.12508862 



-0.4365452849 

0.8996822851 



sin o sin 
cos>^i cos 
cos cos P cos ‘L 
cos 0 

1-1 



- 0,0917015172 

0.8796086600 

- 0,1100290334 

- 0.2017305506 



sxn 



cos y?, cos J3i 



sin^ — + sin*^ 



4 ^] 



V2 , 

sin i 

s1b2i1 



CM)/^ V 

sin2 {^^- A)^/2 
cos Ij, cos B-l S^2 1/2 



cos 






cos 



'sin 0/2 

9/2 

sin ^ 
cos 0 



R sin2^ L ^ gin^ 
6 = 



48°35'34".573 

0.75002937 

0.5625440559 

-38°00'34".710 

-19°00»17".355 

- 0.32564771 

0.1060464310 

0.4948186232 



0.6008650542 

0.7751548582 

50°49'08".668 

101°38'17"o336 

0.9794410574 

-0.2017305506 
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TABLE 2h 



CaJPUTATIQN OF ( A " I-)* ^ 12 ^ 2l« 



(A-l) 



Kl 

K2 

(Kl + 8 ) 
(Ki+ B )Q 

(cVd) 

(C +■ D) sin U 
(F + E) 

(F+E) G 



(Kl H- B ) ^ (C -t- D) sin 0-^ (F-t-E)G 
K2 

237.2388918 
70^19 .5llii5 
238.03280^2666 
U22. 2^123586 
0.028iiii2U9Ul 
0.0278U771 i 65 
0.0l;007it9^31 
0.12875U6570 



Ap« 



( A“ = 

A= L + ( 



U22.U078U83 
70519, 5 Hi ;5 



(.8910182190) ( P ) = 0° 18'20«.866 



A” I*) = 97° Il«09".lli5+0oi8«20'',866 = 97° 29»30''.0H 



cot CX 12 - [tan cos cos A sin A] / sin A 

cot(Xio_ [ (-0.U852216078 ) ( .9776880956)-(-.1303820U28) (0.2100618662 )] 
5.7^Iii63B351 

cotO(i2= -0.UU700709U5/0.991ii638331 = -0.1i5085^6737 
0(l2 = 11U° 16« 06«.607 



cotO( 2 i = [sin cos A “ 003 ^ 2 ,'^®' /^li ^ A 
cot a _ [(-.li365U^28ii9) (-0.1303820U28)-( .8996822851) ( .21ii855706l)] 

0.9911*638331 

cotOC 2 i= -0.1363882066/0.9911638331= - 0 . 137 $^ 8 l 282 
a 21= 277° 19" 56" .501 



cos 



Q 



Qo 

J 3 _ 



TABLE 25 

COMHITATICN OF TRUE SPHERICAL DISTANCE, 6*o* 

(USE A VICE L) 

CCMTOTATION OF GEODESIC 

= sin J^/ sin + cos /^i cos cos A 
= (-=.0917015172) + (,8796086600)(-.1303820it28) 
= (-.0917015172) -(.11U68517UO) 

= -.2063866912 
= 101° 5U* 38 

= 1.778677025 radians 



cos 



cos Uf cos 
sin A 
sin Qa 
cosp^ 



^COS^I cos ^2 a] /sin ft 



0.8796086600 
0.9911t638331 
O. 978 U 70508 O 
0.8912891769 



maximum latitude 
of geodesic 



Ao 



= 1 + (e 



VU) sm^J3o 



- ( ) 
^ -5T ^ 



sin^ 



% 



e 2 

sin'^A. 

3e‘'U/61i 

sinU 

Ao 



O0OOI6920U26 

0.2056036030 

0.0000021Ii73 

O.Oli22728Ul6 

1.0003it77993 



sip2 p, 
e'\l6 
sinli /3, 

R/ o 



go = (® sin2^^ _(e'' Vl6) sin^^o 

O0OOI6920U26 
0,2056036030 
0.0000028630 
O.Oli22728Ul6 



Cc 

e V128. 

sinUp^ 



- I?5 



sin' 



up 



0.00031477691 



0.0000003579 

0,0U227281a6 

0.0000000151 



cos 2 



(T _ 



sxn Bi s±n.0^ - COS 0 

sin2y9^ 



2 sin P2 


0oU20123732l4 


sin 6 0 


-Ool4365U^28l49 


sxn^ Po 


0.2056036030 


cos ^0 


0.2063866912 


cos^ 2 (Y 


-0,6856357822 


2O: 


133°17'08«.26l 


iiCT 


2660314® 16”. 52 2 


cos hcf 


-0.0598071587 



s 



cos 



TABIE 2$ (Continued) 



= b (Aq Oq sin 2 (/ - C, 



^0 

k 

sin 
cos 

Co 

sin 2 
cos U 



00 
2Cf 



cos 



2 Cf 



sm 



2 0 o cos h Cf 



b 

s 



1.00031^77993 

1.778677025 

1.7792956U71 

O.OOO 3 U 7769 I 

O. 978 U 70508 O 

-0.6856357822 

-0.000233309U 

0.0000000151 

- 0.^0388658114 

-0.0598071587 

o.oooooooooU 

6j356,9ll.9li62 

Ils309,3li2.6228 

6,102.141421 

7,027.5723 



sin 2 0 q 



meters 

nautical miles 
statute miles 



7. CCNCLUSICK 



Upon reviewing the results of the determination of the orbit of 
Explorer I, it can be stated that^ in general, reliable orbit para- 

I 

meters cannot be reduced from only one set of three observations of 
the satellite o 

Based on the results of the calcxxlations of i and 0 , the position 
of the orbital plane in space can b e accurately determined from a few 
observations. To determine the size and shape of the elliptical orbit 
and to accurately orient the orbit in its plane requires a series of 
precisely obseived and reduced coordinates. As Moulton states [U,P<>Uij J 
the difficulties arise in finding a, e, and w. As mentioned previous- 
ly, these three parameters determine the size, shape, and orientation 
of the orbit in its plane. To determine these three parameters satis- 
factorily, the whole time interval of the three observations should 
be divided into two nearly equal parts by the second observation. If 
this is so then Cf ^ and Cf^ will be nearly equal. With the time inter- 
vals equal or nearly so, the ratio of the area of the sector to the 
area of the triangles (see Figure 5) will be more reliable and the 
series which is used to determine the auxiliary parameter p will 
converge more rapidly. Additional terms can be used in the series 
to determine p. 

In an analysis of the solution for ui, ug, and u^ (derived from 
p), upon which (7^, (f 2» snd (f ^ are based, the reader should note the 

value of each u determined by its sine and cosines 

U]_ U2 

sine value 21*8° 13 » 02" .933 2^3° 13»50«.317 2^^° 3U’’25”.033 

cosine value 2U8Q 13 » 02". 7^6 2g3 13»50”.30g 2g5o U6»17”.8y4 

difference 0".177 0''.012 ll>g2".787i! 



76 



There is a close agreement between the value of the angle determined 
by the sine and that detennined by the cosine for and U2 but there 
is an ll»52”»78l disagreement between the value for U3 derived from 
its sine and cosine » The value used for U3 in the ensuing calcula- 
tions was an average of these two which is 2^5° i|0»21''<,U2l4. This could 
be the source of the error existing between the computed a, and w 
and those calculated from the 82 observations o I believe that this 
can be attributed to the uneven time interval of the observations s 




22^9“lOS.UiO ) 

22 31 OU, U20 ) 

22 32 03 .100 ) 



l“5il®.280 

0“58s„680 



and the unrealistic observed declinations s 



Obso #1 Obso #2 ObSo #3 

U 70 22« S 51 ° 06« S 1*9° 25> S 



It seems logical to assume that the declination should either increase 
or decrease progressively during the short intervals of the observations 
particularly during the same transit of the satellite o There is no 
such progression in this caseo 

The assiimption that the orbit is correct is another source of 
error. The accuracy of the orbital elements predicted for the time 
of the observations firom the unknown statical depend on the accuracy of 
the Initial values of the parameters determined from the 82 observations 
The standard error of one observation of 5’ »U3 undoubtedly lessens 
the reliability of these values 5 especially the mean anomaly, M, which 
is probably the crux of the entire problem. 

The assumption that the observed coordinates of the satellite 
and the computed orbital elements are referred to the same mean coordi- 
nate system also influences the result. As stated in 9aPo2 on the 
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orbital elements $ 

the reference plane is defined as the true equator of 
the dateo The origin of right ascension is a line 
shifted from the mean equinox of the date by an amount 
equal to the precession in right ascension between 
19^0 oO and the dateo 

Given below is a formula with which values can be obtained to correct 

> 

the right ascension given in the orbital elements, in a right ascension 
referring to the mean equinox of the dates 

(T)= Q° (Doi) + 3°..508 X 10-5 (Smithsonian Day - 33286) 
where DOI indicates the values determined by the Differential Orbit 
Improvement Programo This information was not available. 

It is apparent that the methods and procedures demonstrated 
herein to effect the tie between the two points on the surface of the 
earth can, with refinements, be employed in practical applications® 

The success of this applicatioi depends substantially on the charac- 
teristics of the satellite itself. In fact, the only satellite that 
would be practical for precise results would be a geodetic satellite - 
one that would have the foUcwing characteristics j 

a» Have a nearly circular orbit or one that has a 

very small ellipticity - this would permit reliable 
predictions to be made in its position and eventually 
an ephemeris could be made available | 

b» Be spherical in shape - so that its motion in space 
would not affect its presentation svirfacej 

c. Have a perigee of approximately 1,000 miles to 
minimize the effect of atmospheric dragj 

do Have a large inclination angle so that it would 
be visible from most latitudes j 

e® Have a flashing beacon that would be visible at 
its long range® 
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